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CEAPTER - I 


GEIffiEAL I1T2K)IXJGII0N 

1.1 IM'SQDIJCTIQIv 

The term "Bio-flmd D^chanxcs" refers to that study of 
biological systems where the concepts of fluid mechanics, relevant 
to physiological flows such as blood flow through the cardiovascular 
system, flow throu^ reproductive organs, air transport through the 
bronchial tree etc., are used. Since the beginning of the civili- 
zation, man has been trying to unveil the mysteries of the fumtioning 
of the human body. In doing so he has used various theories and 
laws contributed by our great scientists such as Harvey, ikiler, 

Hewton, Poiseuille, Borelli, Helmholtz, Boyle, Hooke, Young, Pick 
etc, (Pung : Appl. Mech. Review (i960)). 

Since the jiiysiologxoal systems are very complex, studying 
them Idealistically is a very difficult ta^. Several simplifying 
assumptions have been made in the past and even now, to study them 
analytically. That is why in certain areas of physiology the 
experimental research is much advanced than the corresponding 
analytical investigation. 

In this thesis an atten^jt has been made to explain the 
mechanisms and causes of functioning of the following systems by 
making a mathenatical model which takes into aeoaunt the behaviozr 
of the blood such as its viscosity variation, non-Hewtonian character etc. 
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(i) Blood flow throu^ elastic arteries. 

(li) Swinining of spermatozoa in cervacal canal. 

(iii) Dispersion of a solute (drug) in a fluid with variable viscosity 
aM wi 1 ±i non-Newtonian behaviour, 

A brief introduction and survey of these topics are given below 
so that the work presented in this thesis could be seen in its proper 
perspective. 

1.2 BLOCS) ILQy 

Harvey (1615) is credited with the discovery of blood 
circulation which is effected by the pumping of the heart, containing 
four chambers, two on the nght and two on the left known as auricles 
and ventricles. Blood is pumped to various parts of the body through 
the arteries from the left half and brought back to the n^t half 
of the heart via veins. This is called ‘systemic circulation’ . The 
second circulation is called the ‘pulmcnary circulation* in which 
the impure blood ireceived from vanois parts of the body by the 
n^t portion of the heart is sent to lungs for purification, From 
lungs it goes to the left half of the heart for being pumped to 
various parts of the body again. 

Aorta, the main blood transporting tube emanating from 
the left part of -the heart branches out at seveiral points fecilitating 
blood supply to various parts of the body. "All the major blood vessels 
have the property of longi-tudinal decrease in the area of cross-section 
known as ‘tapenng*. Towards the e»i of the arterial system, we have 
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the smallest axteiy called ’artenole’ which again subdivides into 
large mmher of capillaries. In capillary region, the exchange of 
osygen and carbon dioxide between the tissues and blood takes place. 
Capillaries 30 in together to become ’venules’ which in turn join to 
become still larger vessels ’veins’ , The veins from the lower portion 
of the body unite to form ’inferior vena cava’ and those from the 
upper portion of the body ^oin together to f oma ’ superior vena cava* . 
The general size of all the blood vessels and their wall structure 
are described in Burton (l965)f Li^tfoot (l974)« 

Blood is the transpoirting fluid which carries with it oxygen 
and various other mtrients for the body. It is a suspensicai of red 
cells, white cells and platelets in plasma. There are about 5 x 1C^ 
cells in a nnlli-litre of blood. About 5?^ of these are platelets and 
l/65^ are white cells. Middleman ( 1972 ) . 

H-asma is a fluid in which cells are suspended to form the 
whole blood. It consists of 9C^ of water and 79^ of proteins, 

Middleman (1972). The major constituents of plasma proteins are 
Albumin (4.59^ by wei^t), Globulins (2,^) and Pihcino^n (0,59^), 
Lightfoot ( 1974 ). It is established and generally accepted that 
plasma is a Hewtonian fluid. Careful tests have been conducted to 

-1 

determine its viscosity over a wide arange of shear rates, 0,1-1200 s 
and it is found that it behaves like a Newtonian fluid, Copley and 
Scott-Blair (i960). Charm and Kurland (1962), Cokelet et. al. ( 1965 ) 
and Merrill et, al. (l964» 19^5) • There are, however, some reports 
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regarding the non-Newt oru. an hehavioiir of plaaaa hut McDonald (1974) 
attributes this to the presence of denatured protein layer at the 
interface of plasma and air xn the viscometer. 

Red cell is a flexihle hi -concave disc having a radius of 
about 4»5v and thickness 2,5y at the periphery and I.Oy at the centre, 
Burton (l965)» Lightfoot (l974) and Middleman (1972). Its mam 
function is to carry 03cygen from lungs to tissues and get carbon 
dioxide from tissues to the lungs. The volume percentage of red 
cells present in the blood is called •haemotocnt' . For a normal 
healthy individual the haemotocnt value willbe about 4^» It is 
found that an increase in haemotocnt causes an increase in the 
viscosity of blood. At low ^ear rates red cells form a^regates. 

A long chain of red cells thus forming an aggregate is called a 

roulex. Fibrinogen, a plasma protein is responsible for forming 
roulex of red cells. The actual mechanism with which rouleau are 
forned is still unknown. The whole blood exhibits non-Newtoraau 

behaviour particularly when it is flowing under low shear rates. 

m 

Several models have been suggested for the non-Newtonian behaviour 
of tl» blood. 

It IS suggested that blood behaves like a powerlaw model 
for low shear rates. Charm and Eurland (1962), Hershey and Cho (1966), 
Huckaba et, al, (1968). But it has been found that this model does 
not hold good for shear rates more than 20 s” , Charm aaxi Kurland ( 1965 ). 
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It IS also olDserved that at low shear rates hlood exhiMts 
•yield stress’, Cokelet et. al. ( 1965 )* fluid capable of resisting 
motion until the shear stress attains a finite value is said to possess 
•yield stress’. Reiner and Scott-Blair (1959) suggested that blood 
behaves like a Casson fluid. It is shown that the Casson equation 
fits well over a wide range of shear rates 1 - 100,000 s , Charm and 
Kurland (1962, 1965) and 0.1 - 20 s \ Meml et. al. (1965» I964)j 
Cokelet et. al. (l963)» It is therefore generally accepted that blood 
behaves like a Casson model fluid, particularly at low shear rates. 

In major arteries which are elastic the blood flow is 
pulsatile. Euler (1775) was the first to consider the blood flow 
through an elastic tube by considering the blood as an incompressible, 
inviscid fluid. A comprehensive treatment of wave propagation -in a 
VISCOUS fluid in an elastic pipe can be found in the worics of 
Womersley (1957» 1958) and McDonald (1974)* Anisotropy and longitu- 
dinal and radial constraints of the wall are considered by Atabek (1968) 

and Itirsky (1967) m analysing the problem. Cox (1968) investigated 
the effect of visco-elasticity of the wall on the flow behaviour. 

Whirlow and Rouleau ( I 965 ) discussed the periodic flow of blood 
through thick-called elastic tube. Attinger (l964)» Middlonan (1972) 
dealt with the pilsatile flow of Mewtonian fluid in an elastic pipe. 
Lightfoot ( 1974 ) discussed the pulsatile flow of a Newtonian fluid 
in an elastic pipe and through a branching circuit. Rodkiewioa et. 
al. ( 1971 * 1975 ) experimentally investigated the flow behaviour at 
a braiaohing point in the case of rigid pipes. Also Ferguson et. al.(l9l2) 
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studied expeirxinentally the effect of hranching in rigid pipes on 
flow behavioiir of a fluid. 

It nay lie concluded that the blood transport in the cardio- 
vascular system, in general, depends upon the folloring, Fung et. 
al. ( 1972 ). 

(i) Nature of the blood, its viscosity, its non-Newtonian 
character, cell concentration and distribution exc. 

(ii) Blood vessel, its geometry, size and shape, curvature, 
branching, tapemng etc.; its elasticity (flexibility), 
porosity etc, 

(iii) Nature of the flow such as pulsatile, laminar, turbulent, 
micro-circulation etc. 

In this thesis the blood flow through elastic arteries 
has been studied by considering its non-Newtonian behavicur, Ihe 
effects of branching have also been investigated. 

1,5 MOIION OF SFEBMATCgOA 

Spermatozoa are living cells produced in male reproductive 

system. Each spermatozoon has a head and a tail. Its length is 5 2 : 10 ^ 

-5 

cm and tail diameter is 10 cm. It can propel itself in forward 
direction by sending waves of lateral displacement down its tail. 

Its speed is 50yseo, Smelser et, al, (l974)« 

During coitus sperm cells are deposited in the vagina, They 
Efiist move through the bbicus filling the cervical canal to 3»ach the 
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o-viduct. As soon as ovum xs released from the ovary, it is sucked 
into oviduct and it starts moving towards uterus. Because of the 
short life span of the ovum, the fertilisation should occur in 
oviduct itself, iTransport of sperm from vagina to oviduct is too 
rapid to be accounted for by the motility of the sperm. The role 
of oervical mucus in helping the sperm to move fast has been a matter 
of active research and controvert. Several mecMnisms for this 
have been suggested and the woik: done in this regard has been 
reviewed, Sobrero (l963), Davajan et. al. (l970), Moghissi (l969, 1971 ), 
Hoghissi and Blandau (l972), Smelser (l972). 

Odeblad (l962) proposed a transport mechanism of the ceryical 
mucus* Ife studied the physical structure of cervical mucus and found 
it to be a suspension of maeromolecules in a water-like liquid (l959, 
1968). Its viscosity is 0,03P, Odeblad (l962). At mid-cycle period 
these molecules group together ard align themselves near the wall of 
the cervical canal. In the luteal phase, the mucus is like a close 
mesh, Shis structure of mucus in mid-cycle period is supported by 
Elstein et. al, (l97l), Bavajan et. al. (l97l). 

Odeblad (l962) suggested that the long chain of molecules which 
aligned themselves along the canal, vibrate fecilitating quick transport 
of the ^eim, 

laylor (l95l) was 1iie first to consider the motLoa of spermatozoa 
mathematically, asdnniing of micro -or^nisms in an unbounded meditaa 
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was also discussed by Hancock (l953), Gray and Hancock (1955), 

Srokaw (l970) A sin all amplitude solution of th,e swimming of a 
sperm in a bounded medium v/as furnished by ‘'Reynolds (l965) A 
long wave length solution was discussed by Slack and Lardner (l974) 
Ctonsidering Odeblad's liieoiy, Haelser et al (l974) investigated 
the s w imm i ng of spem in an active channel in viiich the walls of 
the channel also did vibrate Ii^ thill (l975) presented a general 
theory of swinming of micro-organisms at low Reynolds numbers while 
pironneau (l975) discussed the motion of uniflagellated micro- 
organisms using resistive force iheoiy 

Ghe actual mechanism of the dynamical interaction of Ihe sperm 
and the cervical canal is not well understood Here an attempt is 
made to explain -this by using the concepts of lubrication theory, 

Cameron (l962) 

1 4 HISPIRSIOH 

Ohe dispersion of soluble matter in laminar flow has biological 
applications such as drug distribution in the body Ihe diffusior/ 
dispersion process in fluids has been studied by many investigators 
under steady and unsteady conditions, Crarfc (l975), Craiic and lark (l968), 
Crank and Gupta (l972), Gumming et al (l966), Gill et al (l970,7l) 

One of Ihe simplest meliaoas to study ihe dispersion was suggested 
by Sir Geoffr^r iCaylor (l953) viio investigated the dispersion of 
soluble matter m a solvent flowing under laminar conditions m a 
circular tube He found 1iiat relative to a plane moving with the mean 
speed of the flow iiie solute diffuses with an apparent diffusion 
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coefficient R^^/48D #iere R is the radius of liie pipe, TJ is ihe mean 
Speed and D is Ihe molecular diffusion coefficient Taylor imposed 
certain conditions vAiich v/ere later raaoved by Aris (l956) 

Numerous authors utilised the Taylor's approach to discuss the 
dispersion in Newtonian and norwNewtonian fluids by considering different 
situations Pan and Hwang (1965) considered powerlaw fluid vdiile Ran 
and Wang (l966) investigated the diffusion in Bin^am and Ellis models 
ReineJvPhillippoff model was studied by Oioshal (l97l) Shah and Cox 
(l974) solved the problaa for E^yrirg model using Taylor's approach 

It may be also noted that, in a wide range of problems, liie solute 
may have homogeneous and heterogeneous reactions wi1h the solvent during 
the process of dispersion. Eor steady state eondi+ions, some authors 
Katz ( 1959 ), Walker (l96l), Solomon and Hudson (l967) considered iiie 
effects of homogeneous and heterogeneous reactions on dispersion 
Recently Gupta et al (l972) discussed the effect of homogeneous and 
heterogeneous reactions on liie dispersion in the laminar flow of a 
Newtonian fluid, between two plates 

In all these investigations, tie effect of viscosity variation 
across +he tube has not been taken into consideration. However, for 
fluids of suspensions such as blood, the viscosity varies across the 
tube, Haynes et al (l959), H^nes (i960), Wiitmore (l968). Middleman 
(1972) Also certain biological fluids such as blood, behave like norv- 
Newtonian fluids under certain conditions, Cliaim and Kurland ( 1962,1965), 
Cokelet et al (l963), Merrill et al (1963,1965), Wiitmore (l968) 

Keeping the above in view, an attempt has been made in this 
thesis to study the effects of viscosity variation and 1iie non-Newtonian 
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character of txie fluid on the dispersion by taking into account the 
homogeneous reaction of "the solute wxtn the solvent 

1 5 SUMM.flRY 

Bais tne^is consists of siy cnapters Chapter I gives the 
general introduction and a survey of 1he work related to blood flow 
through elastic arteries, motion of spermatozoa through the cervical 
canal and liaylor diffusion "through noruFewtonian fluids 

In Chaptex II we have investigated the effect of elasticity 
and brandling on flow characteristics of the blood by taking into 
account its non-llevstonian behaviour It has been pointed out that 
the behaviour of the blood whxdi is non-Fewtoniaa in smaller vessels, 
together with frequent branchings in the pre-arterxolar sections 
causes greater resistance to flow and are mainly responsible for a 
major pressure drop in this region 

In Chapter HI, a generalised Reynolds equation for the motion 
of epematozoa in a channel in a two layer fluid with interaction of 
the wall is derived Die interaction of the wall has been assumed to 
cause a force proportional to -file shear-stress at "the wall Bie 
following particular cases have been discussed in detail and the long 
wave length solution is obtained in each case, Slack and lardner (l974) 

(i) Motion of spermatozoa in a channel in the presence of 
wall interaction 

(ii) Swimming of spermatozoa in a two layer fluid between 
parallel plates 
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It IS fomd that the propulsive velocity of the speim 
increases in the above cases 

Chapter IV deals with the feylor's dispersion of a solute 
with homogeneous and heterogeneous reactions in a fluid with 
varying viscosity flov/ing under laminar conditions in a channel 
or a tube Hie effects of viscosity variation on equivalent 
dispersion coefficient have been studied and it has been noted 
that It increases as the viscosity of the fluid decreases from 
the centre of the duct towards the wall 

Ey using iEaylOr’s approach in the last two chapters of the 
thesis the dispersion of matter through non-Ifewtonian fluids has 
been investi^ted by taking into consideration the effects of 
chemical reaction in the fluid The cases of Gasson and Bic^am 
fluid models have been considered It has oeen noted that the 
effective dispersion coefficient decreases due to non-Fewtonian 
behaviour 
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CHiOPiEER II 


IDJff-lIErTOEEAlT BMiAYIOJR OP BLOOD OIHEDTJGh ELASTIC 

JffiTMIES 


2 1 CTTROIIJO!EEOH 

It IS well known that the walls of blood ressels are elastic 
and distensible Elastin fibre which is present in the tissues of 
the walls is mainly responsible for thear elastic natiire To under- 
stand the realistic behaviour of the blood f-^ow in arteries, this 
effect should be taken into account Euler (177*" ) was the first 
to study the wave transmission in an elastic pipe by considering the 
blood as an incompressible inviscid fluid Later FamerslqT- (l957) 
discussed a theoay of pulse transmission and oscillatoiy flow in 
mammalian elastic arteries tarsiy ( 1967 ) s+udied the effects of 
arasotrofy by taking into consideration the longitudinal tethering 
of the wall while Cox (1968) considered the effects of viscoelasticity 
Ube pulsatile flow in an elastic tube has been discussed by Middleman 
(1972), LLghtfoot (l974), McDonald (l974), Ihe effects of ba^anching 
also have been taken into consideration, la^tfoot (l974) 

In these studies the blood has been considered as a Fewtonian 
fluid Howevei' the blood is apart iculate suspension of red cells in 
plasma and it behaves like a non-Hewtonian fluid, Merrill et al (196?), 
Hershey and C3ao (1966), Hackaba et al (l968). Charm and 
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Kurland (l962, 1965), Y/laxtmore (l968), Bergel Cl972) It has been 

pointed out by Ghana and Kurland (l965), Huckaba et al (1968), Hershey 

and Cho (l966) that the blood be characterised as a powerlaw fltiid, 

when the shear rate is less than 20 sec However, Heiner and 

Scott-Blair (1959) suggested that the blood behaves as a Casson laodel 

fluid. It has been shown +hat Gasson model fits well, for shear rates, 

•*1 

1-100,000 sec , Oharm and Kurland (l962, 1965), 0 1-20 sec , Mernll 
et al (1963), for a wide range of haemotocnt values 

It mey be noted that the cardiovascular system is a very complex: 
net work and the behaviour of the blood flow throu^ it depends, in 
general, upon the nature of the blood vessels such as elasticity, 
tapenng, curvature, branchings etc , the non-New+oraan behaviour of 
the blood, the nature of the flow such as pulsatile, laminar, etc 
It IS not mathematically tractable to have a complete realistic model 
with all these complexities and "some degree of approximation and 
restriction of the generality of the model must be adopted to make 
ary progress at all", Bergel (l972) 

dto study the perfomoance of ary vascular bed, the concept of 
peripheral resistance or resistance to flow,viiich may be defined as 
the ratio of the pressure diop to the flux is most useful, Burton (1965) , 
Bergel (l972) In general the peripheral resistance (ki) depends upon 
the geometry of the bed, ibeology of the blood and the flux Q etc 
In a complete vascular bed the total resistance to flow is the algebraic 
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sum of the resistances in arteries, arfeenoles, capillaries, veins , 
connected in series Erom the definition of peripheral resistance 
it can be seen that for a given flax, the ER is proportional to the 
pressure drop Ihe resistance in the proxirnty' of arterioles constitutes 
a major share of the total resistance fhe mean arterial pressure of 
100 mm Hg drops to 30-35 mm Hg in arterioles, FcSonald (l974-) 

Keeping the above in viev/, in this chapter, we investigate the 
effects of elasticity and branching on the peripheral resistance in 
the arteriolar and pre~arteriolar section of the vascular bed by 
characterising the blood as Casson and powerlaw fluids 

2,2 BASIC B3UAfI0ITS 

Consider an axi^nmaetncal flow of a non-Kewtonian fluid in a 
circular elastic artery whose physical oonfi^ration axd the corresponoirg 
coordinate system are shown in the figure (2 1 ) 

Ohe Httvier-Stokes equations governing the flow in this case are 
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She contmuity equation can be written as 
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Id the model m assuae that the displacoaent of the 

wall xs small so that the radial coiiporifc'"t v of the velocity is 
negligible #ien compared to the ayic*l component u, 3calak in Bing (1966) 
Since the Reynolds niarber of the flow is small, [for floras in the 
region of terminal branches Re = 9 and near arterioles it is 
0 03 , IxL^tfoot ( 1974 ) ] ; the inertia terms are neglected, Batchelor 
(1967) 

With these assumptions, the equations of motion reduce to 


§ (2.3) 

o=-a£ (2 4) 

9 r 

Hie shear stress at ary radial distance r is given by, Bird et al, 

(i960) 


^ 2 dz 


(2 5) 


If IS the wall shear, this gives 


R 


_ I ^ 

2 dz 


(2,6) 


On using equations ( 25 ) and (2 6 ) we get 


T r 



Ohe constitutive equation for a 


(2 7) 


non-Rewtonian fluid, in general, can 


be given as 


25 


du 

dr 


= f(t) 


(2 8 ) 


2 ^ 
dr 


{Ohe flvix Q xs defined as 

E E 

Q = 2Tr / r udr = -ir / 

0 o 

Using the equation (2 7 ) and (2 S) v/e get 

irR^ A 2 


(2 9) 


Q=^ / 


3 


fCx) d' 


(2 10 ) 


!Ehe penpheiQl resistance (SR), denoted by X , is defined as 

AP 
^ “ Q 

Tjitoere Ap is the pressure drop across the bed considered 


(2 11 ) 


It can be seen tram, equation (2 1l) that, for a given flux, 1iie 
SR is propoartional to the pressure drop , McDonald (l 974) 


2 3 CASSOR MODEL REDID 


Ihe stress-strain law for Gasson model fluid is given by 


1/2 / du,,l /2 ^ 1/2 

where Tq aoi n are the yield stress aM consistency Hence 


(2 12 ) 


i(-r) [t + T - 2/rr] 


(2 13) 


Using equation (2 13) m (2 10) we get 
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!I3ais eapressxon is also giren in Whitioxe (i968) 
for anall we get 


Appro xxmtirg 



R 



El im inating from, equations ( 26 ) and (2 15) we obtain. 


(2 15) 




(2 16 ) 


where R * the radius at any axial distance z ^ depends upon the 
elastic nature of the arteiy, 

She relation between the pressure and the stretched radius rCz) 
of the tube can be written as, Rashevsiky (l945), Morgan (l952), 

(2 17 ) 

0 

Ytiere E is the Young's modulus of elasticity, R^ is the unstretch ei 
radius of liie tube and h is the wall thicfcneSw 


Differentiating the equatLon (21?) with respect to z we get 


dp 


Bh dR 


dz dz 


( 2 . 18 ) 


■dp 

Eliminating C^) from equations (2 I 6 ) and (2 18), we get 




where 


8 
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a K r *' 2T^ai 


(2.19) 


( 2 . 20 ) 
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b = 2^/2'''^ (2 21) 

Shougb, the equation (2 19) can be integrated it is difficult to 
get R as an explicit function of z To achieve this, we solve this 
equation by iteration with 


Tq = 0 , R = R^(z) 

as starting values, vdaere R^(z) satisfies the following equation 



The initial condition for this equation is 

z = 0, = R3_ (2 23) 

where R^, the radius at the entrance , is given ty [from equation 
(2,17) ) 


R:, = Ro/[i -P:,Ro/ai] 


(2 2 ^) 


Nowf integrating equation (2 22) and using (2,23), we have 


Rl =Ri(l 

where 
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g _ _ 

Ih ttR^ 

Iterating 1iie equation (2 19), we have 


(2 25) 


(2 26) 
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= -[a R 
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(2 27) 
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waere xs the first iterated yalue of P(z) 

Integrating the equation (22?) using the initial condition 


2 = 0 , Ej - 


(2 28 ) 


we have 






For 6 «1 this can be approximated as 




(2 30) 


where 


♦ -fif 

TTR^ 


(2 31) 


If denotes Idle radius at a = a , we haye from eq.uation (2 30) 


“ ®1 t ^ 3 


(2.32) 


ITom eq,uation (2. T7)j we obtain 


^x-Pc=“‘t5^-5^1 


(2 33) 


(where p is the pressure at a =2,)Tidiich, on using equation (2 32) 


can be written as 


Pi ■’ Pq 


(l-<>) R. 


(2 34) 
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!I5ie peripneial resistance, denoted by X, is given by, tsee equation 

(2 11 )] 


X = 


2 it 

Q R, 


[. 


(| *^0 + an/ q/ irR^)^ 


^ "a f 20 / 0 / 77 ?)^ 
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(2 35 ) 


In the case of Fewtoman fluid (t = O) flowing in an elastic pipe, 


the peripheral resistance X^^, for the same flux Q, is given by 
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Ih ir r:: 
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(2 36) 


Ihe corresponding peripheral resistance for a Casson model fluid flowing 
in a rigid pipe is given by (as E-»- ®) 


X* 


2ji r8 
R, 


rSjTZs r~, 3^2 

[j *^(~)+ 2n >^1/ ITR^ 3 


(2 37 ) 


For a Fewtoman fluid flowing in a rigid tube, we get 

&r^ a 




TTE^ 


(2 3S) 


a well known particular case 

Dividing the equation (2 35) by equation (2 36) we get 


L. r ^ 


/ T TT 
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] ( 239 ) 


1 + 


When E-+- op, we get the ratio in the case of ilgid tube as 


30 



(240) 


It can be seen that the nianerator of the tem in ■'•lie square oracket 
of equation (2 39) is greater than the denominator, v;hiGh ingilies 
that 


X > X. 


'N 


Hence the IS is greater in the case of Casson model fluid than the 

Hewtoniau model. Hie same result holds in the case of rigid tube 

also, as seen fcom equation (2 40) itor (fi «1 , we get the IS for 

2 

Casson moael fluid flowing in an elastic arteiy (neglecting ^ and 
higher powers) as 




\ 7 


2 yh , / it ,2 
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(24 1) 


TTEj 


It may be noted here "that as E decreases to and hence X 

increases ■fco Hub ER in the case of elastic -feube is less when 

compared to the oase of rigid 'babe for Ihe same flux. 


Erom "{iie above analysis it m^ be no-ted that Hi increases as 
the elasticily of the wall decreases and it is enhanced further by 
the non-Hewtonian behaviour of the blood Hius, the non-elactac 
na'ture of blood vessels and non-Hewtoman behaviour of the blood are 
responsible for a major pressure drop occuring in the arteriolar 
section of the vascular bed whidi is compantively rigid, 

McDonald (1974), See Eigs ( 22 ) to ( 24 ) 
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In larger Tsssels, walls are more elastic and blood behaves 
like a Uewtonian fluid causing less and less pressure drop 

2 4 EggEOIS OP HlAItGHISrSS IH 3BB YASCTJLAR BED 

In this section we study the effects of non-elastic branchings 

in an elastic tube by considering the blood as a Casson model fluid 

The physical configuration, which is more or less close to the natural 

situation in the arteriolar aid pre-artenolar section of the vascular 

bed, IS shown in Eig (2 5) Here we assume that the main artesrv of 

length 6fi.(0 < g < l) diverges into M bratEhes at the point B Ihe 

main artery is assumed to be elastic and the branches to be rigid, 

as the distensibility of artenoles is small #ien compared to that 

of the mam tube, Burton (l965 ),Lightfoot (l974) Kie radius of Idae 

m a i n tube is denoted E(z) and iiiat of jth branch by Hie flux 

in the main tube is denoted ty Q and the flux in the jlh branch by 

Q., It IS further assume taat all branches are of 1±ie same length 
J 

(I-B)jI aM the pressure at the end of each branch is p^ We also 
neglect bending, tapering and entrance effects, Bergel (l972), 

Krovetz (l965)» 

If denotes the pressure at the inlet and p^ denotes the 
pressure at the branching point, then the pressure drop in the main 
tube can be written freon equation (2 41 ), as follows 

Pi -Pb 


(2 42 ) 
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where is the elongated radms of the main tube Dae pressure 

drop in the gth branch can be given by 




(2 43) 


Erom equations (2.42) and (2 43) we have the expression for the 
fluxes as 


Q = 


J!3i[ISSZs. + -Jl /S. (fiZS.)'/^] (2 44) 

^^9 * St 7 H:, ^ St ^ ^ ' 


a (2,5) 

®3 Ug ^ n-ejt T E^ 'd-sH' ‘ '■ ' 

Using the condition of continuity, we have 

Q = I Q (2 

3=1 

Substituting 1iie expressions for Q and Q fiom equations (2 44) and 

V 

(2.45) in equation (2 46), we can obtain the pressure p^ at the 
branching point from the following equation 

Pjj[1-B+6l^c^3-4.9 e(l-3) (1 " I “3 ) 
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(2 47) 


a » E /R . 
3 3 i 


where 
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aibstituting the value of in equation (2 42), the pressure drop 
in'*the vascular bed is given by the eqviation 

-n -r f ^ 

_ G 128 ^o^ 1 

M ^ 49 Rj^ M 


I3 1 o: 


I a 


1 


d 


. [/i. „V2j (2 45) 

7 R M A S’ 1 - 3 ^ 4 “3 i 


" I ^ 

P(j_“P -1/2 

■vriiich can be easily solved as a quadratic in ( — — ) ' 

Bo see the effects of various parameters, it is further assumed 
that all the branches are of equal radius i.e , 


a- q, = 


= a 


■R 


(2 50 ) 


and liien 


= QAf 


(2 51 ) 


In such a case, using (2 50) and (2 5l) in the equation (2,49) and 
solving for P^“Pq» we get 


Pi "Po 


2n fJLf 

\ “ a^/TT 1®: 


^ [| /^r+2n 


(2 52) 


We notice that the above equation can also be obtained by adding the 
equations (2 42) and (2 43) and using (2 50) and (2 51 ) 
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It be noted here that M and a are connected the relation 


d 



(2 53) 


v&ere d xs the ratio of total area of cross section of '&e branches 
and the area of cross section of the main tube Ihe average value 
of d IS 1 28, McDonald (l974) and it ranges over 0 75» 1.02, 1 29 
in different parts of the body, Caro et al (l97l) 


Erom equations (2.52) and (2 53 )» the peripheral resistance in 

this case, can be written as 
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mi 
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(2 54) 


Wien the mam tube is considered to be rigid, then the corresponding 
El IS given by 



2(l-g)il- r 8 

^ •• rr 



+ 1^ [ I / C’C^Q) + 2n / Cl/ ITE^] 2 (2 55) 

0 

where ot* = E./R and M, ^rtnch is the number of branches* can take 
only integer values. Its minmum value is 2. It can easily be seen 
that, for the range of d stated above. 


d 


> 1 


(2 56) 
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Biuation, (2 56) and the fact a < 1 enable us to concliJde that 

V X (2 57) 

Hence the BR of the irascalar bed increases in the presence of 
branchings 

Erom equation (2 54) it may be further noted that the ER 
increases as the distensibiliiy of the mam tube decreases and the 
number of branches increases and it is eciianoed by the non-Newtonxan 
benaviour of the blood (See Pig ( 26 )) 


It IS interesting to compare the pressure gradients in the 
mam and branch portions of the bed Binding the equation (2 43) 
by equation (2 42), using (2 50), (2 51 ) and (2 53)» we get the ratio 
of pressure gradients in the case when all the branches are of equal 
radius, as 
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(2 58 ) 


idler e pressure gradients m the mam €uad branched 

portion respeotiTely 

It can be seen that the numerator of the tern on the ngbt hand 
side IS greater than the denomi n ator Hence 
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wfaach shows that the pressure gradient xn the biTanch is greater 
than that xn the maxn tube A similar result for Newtonian fluid 
was obtained by McDonald ( 1 974 ) 

1 

It may be seen further that for 0 » we get 


h-ll2.>± >1 


(2 59) 


from which it be concluded that the pressure drop is more in 
the branch portion than in the m ax n tube 


IHius, it may be pointed out here that nature has provided 
frequent branchings in arteriole section not only to fecilitate blood 
supply to sTl parts of the body but also to help blood slow down 
in the capillary section SUrther, the rigidity of the walls, the 
non-Newtonian character of blood help to create greater resistance 
to flow causing a major pressure drop in the arteriolar section, 
McDonald (l974) 

2 5 POWM LAW HPID MODE!! 

In the previous sections, we have investigated the behaviour 
of blood by considering it as a Casson fluid In the following, we 
stxdy the same characteristics of the blood flow by taking it as a 
powerlaw fluid under the same assumptions applicable to the same 
physical sittaation as described in Fig (2 1 ). 


lEhe stress-«etrain law for power law fluids is 
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T = 


m 


^ jn_i ^ 

dr ' dr 


/ du 

Mr 


< o) 


(2 60) 


where m and n are coosxstenc^ ard. flo?/ behaviour indices 
Comparing the equations (2 S) and (2 60), we get 

fCx) = 

Substituting the abore expression of fCr) in equation (2 IO) and 
usiEg equation (2 6) we get 


0 _ r 

2m dz^ 3+1/n 

Solving for (dp/dz) from (2 6l), we get 


(2 61 ) 


dp -2m_ r 3 ti+l Q T-” 
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(2 62 ) 


itollowing a similar procedure as in the case of Casson model fluid, 
we get the dijCferential equation for R as 


dz 


“2ni 


3at1 


Ih E 


3n-1 nir 


q3- 


(2 63 ) 


Integrating the above equation, with the 3Jaitial condition 


z = 0 , R = E 


(2 64) 


We get E as 


K = \ t 1 

where E^^ is given by the equation (2 24) and k is given by 


(2 65) 


6mn£ r q-iU 

. -,3n nir 


Eh B 


k = 


(2 66) 
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!Ehe pressii3?e drop across the artery is given by 


Hi r 

Pi -Po =57 


1 


(1-fe) 

Hence the El is ootained as 


l/3n 


- 1 ] 


(2-67) 


X = 


Ha 

Q \ 


(1^) 


l/3n 


- 1 ] 


(2 68 ) 


Fo r k « 1 , we get the IR for a powerlaw fluid f lo\ ing in an elastic 
tube, neglecting k^ and higher powers in equation (2 68), as 


X =. r 5njl Qj a (2 69) 

<3 ™ 

With the same approxnaation the IE for a rigid tube is given by 


X* = 


2m S, 


Q 




(3]“ 


nu 


Dividing equation (2 69) by equation (2 70) we obtain 


(2 70) 


(271) 

X* \ 

Since Rq R^j it may be concluded from equation (2 71 ) that the IS 
of powerlaw fluid in an elastic pipe is less than that in a ngid 
pipe 'Thi s in^ilies that the El increases as the artery becomes less 
flexible and iixis increase due to non-flexibilily is enhanced further 
by pseudoplastic nature of the blood (n < l) 

2 6 EFFECTS OF BRANCHIRg 

Eor the same geometry as in Fig- (2 5) when the blood 
IS considered as powerlaw fluid, we have from 1iie eouatLon (2 33) and 
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(2 65 ), the pressure drop in the mam arteiy as follows 


Hi r 1 , 

where k is given liy (2 66 ) 


This can be approximated for kg « 1 as 


(2 72) 


Pi 


2m r 3n+1 Q,n 
_3nt1 '■ rnr 


(2 73) 


The pressure drop in the jth branch is given by 
2m(l--8))l 3n + 1 


^b ~ ^0 


R' 


3n+i 


[ 


Q ] 

mr 


n 


(2 74) 


Following the same procedure as in the case of Casson model, we get 


the pressure pj^ at the junction as 

! „-A,n 




(2 75) 


where 


“3 = ® A 

SubstLtutxng -&e value of pj^ from equation (2 75 ) in equation (2 73 ) 


the pressure drop for the system is obtained as 


P. “ = 


R: 


3£M-1 '■ nir *• 


1-3 

3 


] 


(2 76 ) 


The "SB IB given by 



45 


X 


B 


2m£ l- Sw-l -.n 
,3n+1 ^ n-tr ■' 




1 - 3 


a 


^3+l/n,ii 

D 


] 


(2 77) 


If all the brajaches aare of equal radius, tne eauation (2 77 ) can be 


written after using equation (2 53) as 


X 


B 


_ 2 m 2 __ c 30+1 Q-iO 
qj^3o+1^ o. 


[ B + 


Qt1 

IlLIIz^I 

3 0+1 


(2 78 ) 


In case the mam tube is also rigid, we have 


^B 


2a& 


Q 

O 


n+l 

(a*) " 


where d* = Mcd<- 


(2 79) 


We notice from the equation (2 78) that HI increases as the main 
vessel becomes less flexible It is also seen that the HI increases 
as the number of branches increases 


2 7 CXIHCIPSIOHS 

(The effects of branching and elasticity on the peripheral 
resistance for blood flow have been discussed by charactensing the 
blood as Casson and powerlaw fluids It is shown that the elasticity 
decreases the resistance iidiereas the non-Uewtoman character and 
branching enhance it 

It IS suggested that in the prearteraolar section, the 
rigidity of thie tubes, the non-Kewtonian character of the blood and 
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the frequent branchitjgs, increase the resistance to flow, causiog 
a major pressure drop !I3me, the analysis presented here pirovides a 
liteoretical explanation of the ptysiological observation, pointed 
out by McDonald (l974) 



lOMniCIATORE 


Eauio of the total ^ea of cross section of all branches 

to the area of cross section of the main tube 

Young’s modulus of elasticity 

Wall thickness of the blood vessel 

Length of tne bed 

Himber of branches 

Consistency index of powerlaw fluid 

flow behaviour index of powerlaw fluid 

pressure 

pressure at the branching point 
pressure at the entrance of the bed 
pressure at the end point of the bed 
flux in the main tube 
flux in 3 th branch 

Eadius of the tube at axial distance z 

Elongated Eadius at the entrance of the bed (z = o) 

Eadius of the 3 th branch 

Unstretdaed Eadius of the main tube 

Starting value for iteration of E 

first iterated value of E 


axial velocity 
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a 

“j 

3 

Y 

T1 

T 

■^0 

X 

% 

X* 

x| 

AP 


Halnx) of Radius of branch tube to main pipe ^dien all 
the branches are of equal radius 

proporta-onal length of the main pipe to the total length 
of the bed 

^o^i 

Viscosity of the blood in tbe case of Casson model fluid 
Shear st3?ess at ary radial distance 
yield stress of Casson model fluid 
Wall shear 

Peripheral resistance, in the elastic pipe, of a non- 
Eewtoman fluid 

Peripheral resistance of Newtonian fluid flowing in an 
elastic tube 

lenpneral resistance in the rigid tube of a Hbn-Hewtonaan 
fluid 

Peripheral resistance in the rigid tube of a Newtonian fluid 
Periphecpal resistance due to a Ifiin-Sewtonian fluid flowing 
in a branched circuit udiere the main tube is elastic 
Penjiieral resistance of a Non-Newtonian fluid flowing 
in a branched circuit having main pipe rigid 


Pressure drop 
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CHAPTER in 


SV.'niMPG 0^ SPEHI 1 TO 7 OA IN A CEPnCAE CAETAl 
3 1 lUiffiOEUOTIOE 

Swimming of speimatozoa xn human genxtal organ has received 

much attention in recent times, Reynolds (1965) , Shack and Eardner 

( 1974 ' j anelser et al (l974 ) , la^thxll (1974 , 1975 Pxronneau 

(1975) Ihe spermatozoa are living cells (length 5 x 10 cm, diameter 
“5 \ 

10 caaj having their own motility and move towards the oviduct thro-ugh 
the mucus filling liie cervical canal by sending waves ocf lateral 
displacement down the tail or flagellum A malhematical analysis of 
swutELing of spermatozoa was first carried out by HJaylor (l95l) later, 
Hancock (1953) studied the self propulsion of microscopic organisms 
through liquids Reynolds (l965) discussed a two dimensional model 
hy assuming that Ihe amplitude of 1he tail motion is anall as compared 
to the distance between 1he walls Shack acd lardner (l974) has also 
studied this problem uMer a long wave laagth approjamation However, 
none of these studies has tefeen into account the special structure 
of the mu'*U"' and the dynaaical interaction of the cervical canal 

Cervical mucus is a suspension of macromolecules, having higher 
molecular wei^t, in a water ^ikeliquid. The viscosity of the 
suspension is 0 03 P, Odeblad (l959 , 1962) In the luteal phase 

the mucus resembles like a close mesh, with macromolecules spreading 
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all over the canal, having a spacing of 0 03 ym in *nem, Odeblad (1968) , 
ELstein (l97l ) , Ifevagan e+ al (l97l ) In tre mid c/clc period, the 
molecules group togetncr to form micelles 'diich align tnaaselvcs along 
the walls of the cervical canai, Odeblad (l968 ) 

Sie role played by cerviy in spermatozoa transport nas been 
revdesfed by Mo^-.ssi (l969, 1971) It ha^ been pointed out that, ”the 
transport of sperm through the cervical mucus is too rapid to be 
accounted for by the in vitro swimming speed of approximately 50 y sec 
of sperm in cervical mucus", Snelser et, al (l974) Various suggestions 
have been put forward to explain this rapJ passage of spermatozoa by 
keeping in view the structure of the mucus, tne alignment of the 
micelles along the wall, the possible interactions of the sperm and the 
wave generated by it with micelles, Odeblad (1962 , 1971 ) , Davajan et 
al f"l970 ) A hydrodynamcal analysis has also been presented to explain 
the rapid motion of Ihe speim by assuming 1hat, due to the aoove 
mentioned interaction, the wall of the cervix xinder^es a peristaltic 
wave motion, Saelser et al (1974) 

Here, in this chapter, a differoit approach to explain this rapid 
movement of spermatozoa is suggested by taking into account the above 
mentioned dynamical interaction It may be noted that tiae wave generated 
by the sperm mi^t interact vmth the mice'' les and produce a force in 
tne dirootion of the motion of spermatozoa fSiis force may be caused 
by the flexibility, unevenness and the vibration of the micell ee 
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Purtber, due to the concentration of micelle® near the walls in tie 
midcycle period, the viseositv of the muc s is jreater near t.'>e vialls 
in comoou?ison to t lat at tie centre Ibis crates comparitively more 
viscous resi tance near the wall and causes the sperm to swim only in 
the central region of the cerviy Biis resi tance also slows down 
the motion of fluid near +ne cervical wall which in turn causes tne 
sperm to move faster 

Keeping the above in view, we study here the motion of spermatozoa 
by taking into account 1±ie follo;a,ng aspects 

(i) Eie effec-^ of force, caused due to the dynamical interaction of 
tiie nacelles with the swinming spena, by assuming that it is 
proportional to the shear stress at tie wall 

(ii) Sie effect of a more viscous peripneral layer due to concentration 
and aligmt.nt of tne micelles near the wall in the midcycle phase 

3-2 MAMlIAIPlGil, MiiLYSIS 

Consider tie motion of a spermatozoon in the cefrvical canal wiiose 
jdiyeical configuration is shown in figure (3 1 ) Ihe spermatoejon is 
assumed to be an infinite sheet i(ih.ich propels itself forward by passing 
sinusoidal waves down itseli in the oackward direction 

let Vp be the propelling velocity of -liie daeet in 1iie negative 
x-direction In a framv> moving with the sheet, the form of the sheet 


IS given by 






i^«iiiii«iii 

*' , f0l$l 

' 'Ll'. '■L'"! /' 

lllillHiWIiililii 
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2 = b '^in[k (x - Gt)l= 


(3 l) 


where b is tne amplitude and c is tne relocit/ of tne ware In a frame 
moving with a velocity (c-v)in the positive direction , tne sheet appears 
to be soatiorery and the Huid flow stead Dince t e sheet is performing 
sinasoiddl v/aves, the velocity and pressure distrijutio i are not symmetric 
about y = 0 In addition, as pointed out earlier and shown in fi^ (3 1 ) 
a more viscous peripheral layer exists near the walls in the miacycle 
phase Hence, the entire region between the walls is divided into four 
parts two above the sheet and tvso below it 
+ — 

list p , p be the pressure distrioution above and below the sheet 
and (u'^, v|, w"!), (u^, v^, w^, (u", v", wp and (i^, , Wg) arc the 

velocity components of 1he fluid respectively in the regions I, II, III 
and TV neglecting the inertia effects, as Reynolds number is very 
small (10 ^), the equations of motion in the respective regions are 
given as follows 


(I) 


in the region, 


+ t 

9P r ^ 

+ ]1. [ + 




2 ^ 2 + 

", — ,-n^. ^ 


3X 3y 


„2 + 

a V 


aP 


2 4- 

1 

1 ^ + 


ay 


az 
^2 + 
111 
az^ 


0 


= o 


3u av^ i-wl 


ax 


ay az 


= 0 


(3 2-a) 
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(ll) In the region, ^ ^2 
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2 + 


_2 + 

-2 + 

3P 

+ y 

9 u„ 

t--/ 

9x 



a ^2 

4 — ~ 
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. 2 
3/ 

9z^ 
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^2 + 


^2 + 

2 + 
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9 Vg 
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+ 

t 2 ^ 

4* 


4* :r^ 

9/ 
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9 / 
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dz 


”4* 

^ + 




A 

av^ 





2 

2 

9y 
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ax 

+ 

9z 


V 



(5 2-b) 


(ill) In tne region, -h^ < z equatiomof motion and 

continuity are the same as in (l) There p"^, (u!|^, w^) are 

replaced respectively by p , (u^ , v^ , ) 

(iv) In the region, -h^ < z <-^^t the equations of motion and 
continuity in (ll) hold good where p , (ug, v^j w^) are 
replaced by p,(u2j w^) respectively 


She boiindary conditions are 


(x) 


uj 


rT = 0 


z = h^ 


(ll) 


wj = -bkc cos fcz 




s' 9Z 




f^) 

- e' ^3z ^ 


z = ±h^ 


(3 3-a) 


(3 3 -b) 



± w 
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(ixx) 






(3 j-c) 


4 * ,. 4 - 

fi, „ in 

8z ^2 9 z 

where 0’ is t e constant of proportionalxly and is negative because 

•ftie lorce everted by vJj:i»ting nncelles is assmed to be proportional 
to t-ne shear stress at 1iie wall and acts in Uie direction of motion 

of the speim (boundary conditions (n)) boxmdary conditions (in) 

are smply the contimitj of the velocity and shear stress in x and y 
directions at 1he interface of two fluids inhose viscosities are and 




Hbrudrmensionalising the variables and parameters as below. 


= kx, y = ky, z* = z/h^ 
p± = p± k h|/M^ c 

= u|/c, v±^ ^ ^ ^i/kchg (i = 1,2) 

“ V“’ = V'^2> V = V‘=*'2 ’ “ ' V 

h - h/hg, 0 = khg , s; = B' hg/u^ , b^ = b/h^ 

we get the above eq,uations of motion as follows (after dropping hhe 
i*« imtation for oornreniesice) s 
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(l) In the region, < z < h^ , (h^ = b sin x'* 

2 -h 2 *1" 


ap 

[6^ 



a 

a U.J 

ax 



2 

az 

+ 

ap 

4* 




2 + 

3 v.| 

4* 

2 + 


ay 


9x 


ay 
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] = 0 


au’t 3 t^ 3wt 

■ ^ i 

3x 3y 9 z 


(3 4 a) 


(ll) In the region, h. < z < 1 


+ ^2+ ^2+ 

- r r2 ^ ^ ^ ^^2 

„2®„2^ 2 

3x 9y 3z 


] = 0 


- u 


4- ,.2 4" *f* ,.2 -f“ 

+ [g2!j^^e2_:2 


ay 

+ 


ax 


ay 


az 


au- avj gwj 

— £ + — ^ ^ = 0 

ax ay az 


(3 4-b) 


aLmLlar eq.uationa with superscript may be written for liie 
regions III and IV Using long wave length approxunataon, i e 
e < < 1, and assuming an asymptotic expansion for 1±ie solution , 
as in aiack arfl lardner (1974), we get the eq.ua tions of motion 


involving lowest order terms as follows 
In the region, h^ < z < h.^ 
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(3 5-a) 


(5 5-b) 


(3 5-c) 


(5 5-d) 
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where 



Bie mn-dxm.ensionalxsed boundary conditions are 


(x) 


= 0 


on z = b sin x = I 


(3 6-a) 


w^ = -b cos X 


(ii) U~ = V - 1 ia( 


+ { 
= + a (' 
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!!i 

dZ 


az 


Z - + ^ 


(3 6^}) 


Wg = i w 


(m) 
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1 
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az 

3z 


2 4 


(3 6-c) 


*1 

where a = — ihd is negative 

Solving the equations (3 5-^)> (5 5-b), (5 5-c) and (3 5-d) in 
the r^ons above and below the :Sieet and using the boundary 
conditions (3 6), we get 
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ax 

2 + ^ 
z + ^ z + 

(3 7) 
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2 

— 0 + C- z + 4 

(5 b-a) 

+ 


2 + 

W " C- 0 + D- 

(3 9-b) 


v^ere 

0 p i 

U1-A2CX) + (tT - 1) (1 - + 2 ®)] 

A* : 3 (3 10-a) 

(1 +h^ + a) +(y + a) 


4 = y A* (3 I0~b) 

m» *1 9 !P'” o « p 

+ + 2 a) + (ii - 1 ) (l-h^ + 2a)] 

Q± = (3 10-c) 

(1 + + “) + (v - 1)(1 - + a) 


0* = p- 0* 


+ 1 2 4. 

= - 1 - — - — n - i: h 

1 ‘ 2 ° 


= T^-1-V aF ^2 ^ (1 + “) 


1 si*"" -2 i 

J _ 0* Ii 

o 3 y o 10 




--y 4+0) +0* (1 + a) 

ay '2 2 ^ 


(3 1 O-d) 

(3 11 --a) 

(3 11~fe) 

(3 11-c) 

(3 11 -d) 


It may be noted ttat the equations (3 l 0 -a) to (3 11 -d) contain taie 
unknown velocity V-p, ^ich is to be determined by using the equation 
of force balance on the sheet 
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®ie flux Q , Q ■ in liie x and y directiois aboire and below the sheet 
are gxven oj 


= / u dz + / dz 

h 

1 


b sin X 
-h. 


b sin X 


= J dz + / u. dz 


"x ' 2 

-1 


-h 


1 


1 


= J ^ dz + / ^2 dz 


b sin -X 




b sin x-. 

Qy = / ^ V" dZ + / ""l 

^ -1 -h^ 


(3 12-a) 

(3 12-b) 

(3 13-a) 
(3 13-b) 


Substituting the expresaiom for Telocity components from equations 
(3 7 ) to (3 9 -h)and integrating, we obtain the fluxes as 


^x~~^ 

F 


1 

12 


3 P- 

ix“ 


F" 


(3 14) 


Q± 


1 

12 


9 P 

W 


(3 15) 


where 


F" = (1 + h^ + a) + (iT - 1 ) (1 _ + a) (3 16-e) 

G- = [(1 + (U - 1) (1 -h^)^]^ 

+ 4 (y - 1) (1 ? h^) (1 - h^) (h^ + 

+ 4 ap [(1 + + (ii - 1 ) (1 - h^ J 


(3 16-b) 
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H- = -^ [(l ? + (u - l) (l - ] 

- (l + h^) [(l + + “) + (y - 1 )(l - + a) j (3 16-c > 

Integrating tne continuily eqviation between z - and z = 1, we get 


3 (0+) + f. (Q+) = 0 
3 x ^ ay y 


(3 17) 


Substituting for Q"^ and we get the equation governing the pressure 
X y 

in the region above the sheet as 


r -1 

3x '■ 


,+ ap 


+ aP 


__ £_r_3 

12 ia S' 


+ |- [ ^] = 0 


(3 is) 


Tiiiere G-^, are given equation (3 l6-a) to (3 l6-c) 

Similarly tie equation in the region below the sheet is obtained as 
follows 
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aP 


3x 
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r 3x 


- 3 ...„, 

ay 
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12 P' 


- aP 


9y 


ax 


E~] = 0 


(3 19) 


where p”, g", h” are given by equation (3 16-a) 10 (5 l6-c) 

It IS pointed out here that these equations may be seen as 
the generalized forms of Reynolds lubrication equation applicable to 
spermatozoa motion, Oameron (l966) 
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Particiilar Cases 


( 1 ) When a 5 .^ 0, we get the follovang Reynolds 

equation in t iO dimensions for tne case of a fluid of constant 
Tiscositj in viiich the sperm is smmraing and having the dynamical 
interaction proportional to tne ^ear stress at the wall, 


A_ 


12 i 8x J 
1 




G- 


9P- 


_L 
± 3y 


P, 


1 



-Aj =0 
P7 


(3 20) 


wher“ 


P* = (1 +iio + a) 

G-* = (1 + h^)^ (4a + 1 + h^) 

Y 

^ (1 + - (1 + \ ) \ 


(3 21 -a) 

(3 21 -b) 
{3 21-^) 


In one dimensional case, the above equation reduces to 


[ 



12P^^ 


] 



3x 



(3 22 ) 


where Q * a constant of integration, as can be seen frcaa equation 
of contxmity in this case represents the flux 

(2) When Os= 0 1 e , there is no dsroamical interaction, the 
equation governing the pressure in Ihe canal when the viscosity of 
the mucus varies as a step function, is gLven by 
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12 



3P' 




(5 23) 


where 


= (l + h^) + (v - 1 '' (l - ) (5 21-a) 

G* = [(l ? h^)^ + (il - l) (l -h^)^]^ 

+ 4(v - l) (l ? h^) (l - h^) (h^ T h^f (5 24-b) 

4 = ^ [(1 ?h^)2 + (p - 1) (1 -h^)2] 

- (l + h^) [(l + h^) + (u - 1 ) (l - h^ ) ] (3 24-c) 


In one dimensLonal case 1±iis equation beccmes 


Z±.f2 

„± ^ 3x 

12 Pg 


ap- 


(3 25) 




vfcere Q is again a constant of integration and rejr^ents the flux 


It IS interesting to note Itoat whoi ot= 0 and. V.j = both iiie particular 
oases obtained abofve reduce to 


(1 + Jv ^ > 


(3 26 ) 


which coincides witii the equation obtained bj Shack and lardner (1974) 


In wnat follows, we consider one dimensional case of swinniing of 
apennatozoa 



6d 


Since the saeet is self propelling the resultant fluid 
force acting on it nust be zero !Ehis is expressed as 

/ - 0?“) ds = 0 (3 27) 

***" 

where T , T are the fluid forces on the upper and lower portions 
of the sheet Substituting the expressions for shear, from stress- 
strain l^w and taking the long ware length approxiiaation (& « 1 ), 
the equation (32?) reduces to 




-TT 


[p] 


dx = 0 


IT 3u 

f I [ 1 , + b cos X [pl|dx = 0 

j I'-az z=bsinx 


(5 28) 
(3 29) 


where [ ] denote the difference of the enclosed quantity evaluated 

3u 

above and below tne sheet Appropriate expression focr [ 
be substituted for getting the equation of balance of fluid force 
on the ^eet 

dlie rate of work done by the swiaaang sheet against viscous 
forces per unit thickness is 

where is Idie unit normal to the sheet, u^ the velocities of 
particles on the sheet in the fixed frame, and are the fluid 
stresses on the upper and lower sides of the sheet 


I 
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aibstituting for the shear stresses and taking the long wave 
length approxaiBtion we get, 

TT 

W = / - h [p] cos X dx (5 3i3 

“tt 

3 3 EgPECr Qg -nPUmiOML IlffiERlCglOI 


In this section we Investigate the effect of dynamical interaction 
of the micelles at the walls on the motion of speua in cervical canal 
Eie differential equation for determining the pressure is given hy 
frcan equation (3 22), 


dp=^ = 
dx 


-12Q (a + h^) 

(4a + 1 + h^) (1 + h^)^ 


6 T„ (1 + h„) - 12(0 +1 + h^) 

X O O 

( 4 « + 1 + hQ)(l + h^)^ 


(3 32) 


where the parameter a determines the dynamical interaction and is 
a negative constant pointed out earlier. 


Integratxr^ equation (3 32 ) and expandixg the resultant in 
terms of a(|al « l), we have after neglecting a and its hi^er 
powers, 

2 . N 2 


dP, = - 


- . + 2 ) 3a(3b 4 2.) i 

- (i?p ’ 

+ 12 T], 1 [ 5T75 - g^S/j l 


- 12Tf[ 


2 3 a (b^ + 2) , 




(3 53) 


WWtrtt - VC'^) 
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The force balance equatmn (5 29) xn this case is given by 


I I ] V cosx[p]Viix=0 

■> l’-3z-‘z = b sin X J 


(3 29-a) 


Substituting for [-^3 from the equations (3 7) m equation (3 29-a), 

aZ 

we get, after sxE^jlification, 


8tt 


/ (l+a)^ - 


2Ap + J b sin X [-^1 


d3C 


IT (i-h ) apYa-x; (i+b ) a,p/dx 

f r ^ Q 

“3 a+1-h ^ 0+1+ n 

—IT 0 0 


] ax 


(3 34) 


Substituting for ana evalxiating the integrals, we get fron 

2 

equation (3 34), after expanding in powers of a, and neglecting o and 
higher powers. 


Sir T [ \ Vp 

^ (l-b^)V2 






(3 35) 


Solving the equatiom (3 33) and (3 35) for Q and V^, we get 

6b^ o 

1+2b^ (l-b^)(l+2b^)^ 


(3 36) 
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4Tr 


3b 




a 


(2-b^) 


6 a b' 


477 (1-fa^)''/^ 


(l+2b^)^ 


(3 37) 


When a = 0, there xs ro reaction at the wall 3he result coincides 

with that of Shack and Ijardner (l974) Since a _< 0, in the free 

channel case, when Ap, = 0 , we notice that the effect of dynamical 

A 

interaction is to enhance the propelling velocity of the speim by an 
6a b^ 

amoimt tt-o 

(l+2b^) 

Substituting the expression for [^1 and using the Tsiues of 
Q and Yp, we get the work done, from equation (3 31 ) aa follows 


W 


(1^^) 


2477 b 

172“ 


(l+2b^) 


d-b") 


72a b 

2^ 




(3 38) 


For small amplitudes with aPj^^ = 0> expressions for Yp, Q and W can 
be approximated as, 


2 2 
Y^ = 3b - 6 Ob 

Q = -1 + 3b^ - 6a b^ 

2 2 
W = 2477 b - 720 h 


(3 39) 


Which agree with Shack and lardner (■’974) , Tffcen a = 0 

It IS seen here that since a 5 0, both Yp and W increase as 
the force due to dynamical interaction increases Ihis increase is 
further enhaixjed by the increase in the amplitude of the wave motion 
(i e o) These results can also be seen from the figures (3 2) and 

(5 3) 
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5 4 EgEBGT OF pmaHEEAIi lAHE 


In the pnerxous section m hsve seen the effects of dymimcal 
interaction at the walls on the motion. In this section, we study 
the effects of a more viscous peripheral xayer on the swimaing of the 
sperm, in a channel She differential equation {3 25) for pressure 
distribution in tnis case can be rewritten as 



_2 a£. 
4 ' 



where G-*, are given by equation (3 24--a), (3 24-b) and (3 24-c) 
After substituting the valiie of Fg, Gg, in the above equation, we 
get, 

0 [1 + b sin X + (y - 1 ) (1 - h.^ ) ] 
dx *" D 

6 CTp“2)[(l?b sin x)^+(v-l)(l-h.j)^]-12 (l-h.j )(u-1 )(h^+b sinx) 

+ - 

where 

D = [Cl + b sin x)^ + (W“1) 

+ 4(v-l)(l-h^)(l+b sin x) (h^ + b sin x)^ 

2 

fe assume the pepiphsral layer to be very artall so that (l~h^) and 


higher powers can be neglected In such a case, we get, 
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“12 Q [(l+b sin z) + a^j 

•>* X — 

(l b sxn xf (l + b sxn x 4- 4 ) 


6 (Yp~2) (l-fb sxn x) - 12 ct 

(l-fi Bin z) (4a^ + 1 + b sin x) 


in which is given by 


(3 4l) 


= (ii - l) (1 - h^) (5 42) 

Since y <1, vre notice tMt 

< 0 (3 43) 


Integrating the equation (3 41) between -v and w and expanding the 
resultant in powers of oLj by assuming joj | « 1, we obtain, after 
ne^ecting the small older teanas. 


IPj - -12 Q 1. t^^_^2^5/2 


b\g 3a^(3b^ 




) 


+ 12 Yp » [ 


1 




2a^(b^+2) 


1 


- 12 TI [ 


0 ^ 


3a^(b^+2) 

1^' 


(3.44) 


where APj^ = pCir) - pC~ir) 

3u 

Substituting for [ — from the equations (3 7) and (3 lO-a) in 
32 

the equation (3 29) we get the equation of balance of fluid force 
as 
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Sir y. 


= 2&p, + / to sin z [|5 


0/5 ULiW 

>'(1+0,) - b 

IT (l+b sin x) dp /dx (l-b sin z)dp /dx 

“1 / f a +1 - b sir X **" a + 1 + b sxnz^ ^ (3 45) 

-T! 1 1 

-J- M 

Sub stituting fcQ? and and integrating and expanding in powers 

2 

of a^, neglecting and higher powers, we obtain from eqviation (3 45), 




2 “.( 4 b'^ + tSb'^ + 2) 

~ "^24 Q It [ o ~ 0 7 ffy J 

^ (l_b2)V2 (i-b2)V2 


r 2b 

+ 12 T-n If t 

^ (1-b^)' 


a^(l6b^-i-2) 


24 IT [- 


a^(l3b^-r2) 




(3 46) 


Solving the equations (5»44) and (3 46) for Q and Vp we get, 
AP (l-b^)^/^ ^2 Apj^ (6-9b^) 

n 4* * "" ^ j * i"" ”j" TL" ' ' ■ ' ' 


12 It 


1+2b^ 24(l-b^)^'^^(l+2b^)if 


a^(lS - I5b^) b^ 
2(l-b^)(l+2b^)^ 


APj^ (l-b^)'’^^ 3^2 APx b^(6-7b^W) ^ 

s= 4* r> * [ o ■x/o ' o' 3 


If ^ 8(l-b^)^'^^(l+2b^) 


3a^ (8 - 8b^ H- 3b^)b® 
^ (l-l)^)^(l-i2b^)^ 


(5 47 ) 


(3 48) 


I 
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when = 0 , we notice that 


Q 


I2n . 


1t2b 


AP, 3tj2 




47r 


1+2b'' 


(5 49) 
(5 50 ) 


Which are the same as obtained by Shack and Lardner (l974) 

In a free channel case, = 0 Since _< 0, it may be noted 
that the effect of a more viscous peripheral layer is to increase 
the velocity by an amount 


3a^b^ (S-Sb^+Sb"^) 

^ (l-b^)^(l+2b^)^ 


(3 51) 


For small amplitudes with = 0 » 


we get 


Vp = 3b^ - 12o^b^ 


Q = -1 + 3b^ - 9 


(3 52) 
(3 53) 


It may be further seen here that as the amplitude b increases, 
the velocity ?p increases Ihe woik dor® by the sheet against 
viscous forces is given by equation (3^31) oti integration by 

parts, yields 

W = J b sin X ^ 
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aibstituting for[-^] from eqtiations (5 .41) we get W as a ftoction 
of Q arid Usirg the esj^resaxons for Q and 7^ from equations (3 47) 

and (3 48) and sin^ilifyir^ we get 




(l-h^)’’'^^ (l+Sb^) 


2 4 6 8 
2b - 17b +I7b -8b 

^ (l-b^)^(l+2b^) 


— 56o(-j 1d^ tr [" 


2 4 

2 - 41 ) 'h V 

(l-b^)^'^^(l+i 


■(I4.2b^)' 


(3 54) 


for small amplitude (b « l), with » 0,we get 

f = 24 w b^ - 720^ ■’f (5 55 ) 

It can be noted here that the effect of the peripheral layer 
is to increase Tp and W since ^ 8 See fig C5 4) (5 5) 

3 5 CWCOTSKMS 

A different approach is su^gesfed to account for the quioh 
transport of sperm throu^ cervical mucus by t^ng into consideration 
the ihysical structure of the aucus and the dynamoal interaction of 
the cervical canal Following Shack and Lardner^ the expressions 
for the flux, propulsive velocity and work done are calculated 

It has been pointed out that the propelling velocity and 
the work done by the sperm increase due to the dynamical interaction 
of -the cervix asid due to the presence of the pen^exal layer This 
increase is further enhanced by the increase in the amplitude of the 


wave produced by the sperm 
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lK)HEI»CIA!ineE 


P 

P' 

q’ 


I 

T 


ampll'uude of Ibe wave generated ty ibe movxng spena 
wave Telocity 

h = b am x \4iere x = b sin x is -tbe eiiuatxon of the ware 
o 

in the non-dimensiomlised foun 
half of width of the central layer 
half of Width of the channel 
wave nmber 

pressure above tne sheet 
pressure below the sheet 
flux in »-directlon above the sheet 
flux in y-direction above the sheet 
flux in x-direction below Ibe ^eet 
flux in y-dixeotion below Ibe sheet 
fluid force above 1be sheet 
fluid force below the sheet 


(u^,T^,w^) Telocity vector in ibe central layer above 1be sl^et 

(u^v^,w^ velocity vector in Ibe peripheral layer above the dieet 
2 * 2 2 

(u'*«v"tW*} velocity vectcBT' in the central layer below the ^eet 

\ .j» -j/ 

veLooity vector in the peraphexal layer belcw the sheet 


V 
P 

W 

V 


p-mpftl 1 Irg velocity of the spena 
woifc done by the sheet 
viscosity of tbe central layer 
viscosity of •fee penjheral layer 

pi*eBsure dxffereooe in wave lepgth 
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CHiP3SE 17 


EEPECTS OP Y1SX)S1TI YiELiHOI OF DlSPfflSICf^ IF a PUJID 

4 1 IIIIB3IIKHOH 

In biological ^gystaas sxich as -yascular bed, dispexsion plays 
an important role in the distribution of drug and other nutrients 
One of 1316 pbysical parameteisto be studied is iiie coefficient of 
dispersion Yhich measures tiie rate of distribution. Various attempts 
have been made by several authors to investigate the daaracteristios 
of dispersion in fluid <^nanncal situations vfcidi can be applied to 
biological eystems 

It yms Sir Geoffr^ Ihylor (l953,l954) #io first mggested 

a siiiqplafied approach to study liie problem of dispersion of a soluble 

matter in the laminar fl09? of a leirtonian fluid in circular pipes 

It v*ae found "Hiat, relative to a plane movii® with the mean speed of 

the flow, the solute appears to diffuse with an apparent diffusion 
2 ? 

coefficient ?4iere F is the average speed of me flow, B is me 

molecular diffesaon coefficient and R is me radius of me pipe# senior 
iuqposed certain restrictions on his analysis micfa were later removed 
by Aris (l956) Bimexous aumors utilised Taylor's approach to find 
cOQoent 3 ?ataon profiles in the case of non-Fewtonian fluids Ban and 
Hwang (1965) considered powerlaw fluid while Bmi and Wang (l966) 
investigated, me dispersion In Bin^an plastic and Ellis modal fluids* 
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Bie case of Rexriesr-Pfail: 5 paff model fluid was considered by Ghoshal (l97l) 
aiah and Cox (l974) studied the i^rii^ model Crank (l968,l975), 

Crank et al (l972), Cumming et al (l966), GUI and ^Efcarasubra- 
manian (l970, 1971 ) have also studied in detail the steady and unsteady 
state uiffUsion in various cases 

Ifone of these studies has taken into consideration the homogeneous 
and heterogeneous reactions of Ihe solute with the solvent Recently 
Gupta and Gupta (l972) discussed ''the effects of homogeneous and 
heterogeneous reactions on 1h.e dispersion of a solute in the laminar 
flow between two platei^' 

In all Ihese investi^tions the viscosity cf Ihe solvent is 
taken to be constant. However, for many biological fluids and fluids 
of suspensions such as blood, the viscosity varass across the tube, 
Haynes (i960), H^nes et al, (l959), loh (l^9, 1975), Whitmore (l968). 
Middleman (l972) It is therefore desirable to study the effects of 
viscosity variation on the dispersion of a solute in a flowing medium. 

In this Chapter, we study the effects of viscosity variation 
on the dispersion of a solute in the laminar flow between twi pistes 
and oiroular pipes by taking into comideration the homogeneous and 
heterogeneous reactions of the solvent with the solute 

4 2 hlSPlRSlOM aHRCOGS A EEHID ELOW BECTEES IWO HATES, WIOH 
HCMOGMBOTS REAGTIOH 

Consider the laminar flow of a viscous liquid of varying 
viscosity under a constant pressure gradient between two parallel plates 
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distant 2h apart ahe jiQrsical configuration and the coordinate syatea 
are ^own in the Pig (4 1 ) The viscosity variation is tdfcen as a 
known function 

p=p(y) 0<y<h 

and 18 assumed to be symmetrical about y = 0 i e 

»'(-y) = Ky) (4 1) 

The ftiLly developed velocity profile in this case is given by 

teW y ^ (4 2) 

y ^ 

The average velocity, vdiich is defined as 
1 ^ 

u=-2^/udy (4 3) 

.-h 

is given by 

o 

W© assume that the diffusing matter while dispersing, undergoes a 
first order irreversible ehomoal raction in the fluid under isotheimal 
conditions In such a case, the equaticMi for concentration c is given 

by 

||+uU = D^-Kc (45) 

ay 

where D is the molecular diffusion coefficient* K is the reaction 
rate constant* It mey be noted that, in writing the eijuatlott (4 5) 




v/-' 
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we have aseiaaed litiat 

2 „2„ 

2 2 
9X 3y 

Eelative to a plane movxng with the mean speed u of the flow, 1die 
velocity of the fluid is given by 


= u - u = (- -g) h^ [G-q - f(h)l (4 6) 

Tfdaere 

I Mr 

at (4 8) 

Using equation (4 6), the concentration relative to "Hie moving plane is 
given by 


3o , 3c 

3 t 3x 



ay 


fc 


Kon-dnaenstonalising the variables by 


(4 9) 



where L is a characteristic length along x-axis, the diffusion equation 

(4 9) can be written as 

u 2 

1 3 C ^ X ^ ^ 3 o jr e 


(4 10 ) 
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Pollowxng Kaylor, we row assme 14iat partial equilibnum is oatabli^ed 
in any cross-section of 14ie enamel so that liie diffusion eqtxation 
can be written as 


& . c = H - f(n)J 
3n 


where 


Kh 

B 


H 


(JP)(_ M) 

V Ac/V” j r ' 


BL 


2 ''35 


dC' 


and (-l^) is a constant 
35 


The boundaiy conditions are 


(4 1l) 


(4 12) 
(4 13) 


Is. * 

3n 


0 at ^ 1 


(4 14) 


Solrirg the equation (4 11 ) with Ihe above boundaiy condition, we 
get 

n 


c =s A co^ on 


H snai on 


/ cosh at f(t) dt 


+ 122S12! /” at - ^ 


(4 15) 


\riiere 


A = • - / cosh (a - at) f(t) dt 

a sum o ^ 


(4 16) 
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3]h.e volumetric rate at which the solute is transported across a 
section of the cnannel of unit breadtn is 


Q = / c u dy 
_h ^ 

Using equations (4 6) and (4 15), we obtain 

iff (1^) 

Q.. p(^) 


(4.17) 


U 


(4 18) 


whea?e 


F(a) (f^ - Pg ^3^ 

f(t)dt] / coBhCat>[f(t)-Gj dt 


1 t 

Pg = / ^sihh at [f(t) - GqI t / coda as f(s) dsjjdt 
o o 


and. 


1 t 

Pj = / ^coshQat) [f(t) - Sq] r / f(s) sixii as dsjjdt 


(4 19) 


CkHjparing with Pick's law of diffusion the effective dispersion 
coeffKJierEb It* is given 

if (Mf 

1(a) (4 20) 

It is noted 'tiiat this espression is applicable for any given general 
function IJi(y) However, to study specifically the effects of viscosity 
variancion on dispersion, the following relation is assumed which may 
arise in physiological situations* 
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y (y) = p (5 « 1 ) 

O ' 

or 


(4 21 -a) 


v(y) ~ 


(4 21-b) 


Using the equations (4 8)> (4 21— b) xn equation (4 15) we obtain 
the concentration profile as follows 

2 

c = cosh 001 + M„ (T--i?-n ) 

1 2 5 ^2 

+ 5 cosh ati+ Ug sinh on + Mg ] (4 22) 


where 


= 


^ I£D^ W _ 0^ 


ainh o 




{— )(S) 


2UlS 

o 


(4 23-a) 


(4 23-b) 


and 


U. 


1 2 

BL u sinh a 
o 


r 1 .2 2 cosh a , 

%%) [-3 +“5 5 J 

a a a 




2ir 


'2 2 

Iffi u a 
o 


.5 


(4 23 -c) 


(4 25-d) 


iKie flux Q and the effeotire dispersion coefficient are given 
equations (4 18) and (4 20 ) in vfaach f(o) takes the form 


A - 1- [ fJ- -1- + ° . J-’y + fi(i- ^ -L. „ 

no»6) - o o‘^^45 *■ 2 * ^ 2 4 ^ ? 


2 2 
O 


3 a 

4 eotb a 


36 


2o 


4 

a 


4 coseoh o 

3 


)] 


(4 24) 


W ^ 
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Bo see -fee effects of a analytacallv on dispersion coefficient, 
we expand P(®,®) for small a as follows 



a, 6) = ( 


945 4725 


) + H 


11 2290^ s 

4320 " 945 960^ 


(4 25) 


Ab a 0, 

^Lt.|3Kcx,5) (4 26) 

■wfaich generalises the insult of Wooding (i960) for a fluid of variable 
Tiscosity 

\iien = 0, we get "ttie case of a fluid wilii constant viscosity 
for which we have 


i<“) 



ooth a 


a' 



Kiis agrees with the remit of &ipta et al (l972) Ihr a fixed 5 
we notice that p(®,®) decreases as a increases as observed by Gupta 
et al ( 1972 ) Hiis is due to the presence of honKJgeneoua reaction 
of the solute wife the fluid* But it can be seen from equation 
(4 ^) that as 5 increases F(a,6) itKji'eases a fixed a, idaich inqplies 
that the equivalent coefficient of dispersion increases as the viscosity 
towards the wall decreases 


5to see liie effects of o , in general, B(o»fi), given by 
equation (4 24) has been calculated and tabulated in table (4 1 ) 
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a 


2 0 


4 0 


6»0 


Table 4 1 


6 





^-2 

0 00 

0 15133070 X 

10 



^-2 

0 01 

0 15315713 

10 

0 1 

0 16959501 X 

10“^ 

0 2 

0 18785933 ^ 

10"^ 

0 3 

0 20612364 

10"^ 

0 4 

0 22438796 

ICT^ 

0,00 

0 81988283 x 

10-5 

0 01 

0 82984501 X 

10-5 

0.1 

0 91950461 X 

icr5 

0 2 

0 10191264 x 

icr^ 

0 3 

0,11187481 * 

10-^ 

0 4 

0 12183699 X 

ICT^ 

0 00 

0 46725125 

icr5 

0 01 

0 47296614 

10-5 

0 1 

0 5 2440040 ^ 

i£r5 

0,2 

0 58154956 x 

10-5 

0 3 

0 63869873 

10-5 

0 4 

0 69584790 x 

10-5 
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a 


F(a,S) 


0 00 
0 01 
0 1 

8 0 

0 2 
0 3 
0 4 

0 00 
0.01 
0 1 

10 0 

0 2 
0 3 
0 4 


0 292544^ X 10"^ 
0 29614129 X lO”^ 
0 32850888 x id"^ 
0 56447287 x lo"*^ 
0 40C43686 x 
0 436400® ^ 10“^ 

0 19788888 x 
0 20053066 X icT^ 
0 22230665 x 
0 24672443 ICT^ 
0 27114220 X ICT^ 
0 29555997 x io“^ 
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It can be seen from table (4 1 ) that as ct increases P{a,5) 
decreases for a fixed 6 Ubese results are same as obtained 
analytical!/ 

4 3 EFEDGIS OP CCSIBIHED HCMOSEIEOPS AHP HET5EOSEHKKJS REA-CflOff 

•Considering the effects of chemical reaction taking place in 
the fluid and at the walls the diffusion equation after non-dimensional, 
ising the -variables and following E^lor (l953)» faeccmes 

-2 „ 

- a c =H [G^ - f(n)] (4 27) 

3ti 

where H, G^, f(n) are given by equations (4 15), (4 7) and (4 S) 
respectively 


31he boundary conditions are 

± gc = 0 at Ti= i1 (4 28) 

an 

where 

6 = sh 

and s is 1±i.e heterogeneous reaction rete constant 

Solvu^ equation (4 2?) with boupadaiy conditions(4 28) we obtain 

tl H G 

c w A cOEfo otri - / f(t) sirh a^*t) dt - — (4 29) 

a 0 

where g g 

H [S + + -f)] 

4 a 

A 


asihh o + g coEh a 


(4 30) 
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•} 

= / f(t) cosh (a - at) dt 
o 


and 


1 

So = I fCt) sinh (a at) dt 
0 

{Phe volmetric rate at whicri 1±ie solute xs transported across a 
sectxon of the channel of unit hreadlii is given by ea nation (4 17 ) 
aabstxtutxng for c and from e<iuatxons (’ 29) and (4 6) we get 


Q = 


(|f) 

D 


G(a,e) 


(4 31 ) 


whei*e 


G(a,p) = G^(a,e) + a^(a) 


(4 32) 


S- G 
S + p(-Z + 




1 


2 . £ ) 


^ J 


a sinh a + 6 cosh a o 


[f(t) - G^] cosh at dt 


and. 


SjC-) 


1 t 

= j - oj [ ^ f(s) sinho(s-t)ds]Jdt 

o o 


CcHiq>aring with Mcfc’s law of diffusion, the effective dispersion 
coefficient is given by 


» = 8(a,8) 


(4 33) 


where G(a,8) is given by equation (4 32) 
As in Ihe last section, we assume 
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]i (yj = e""^^ ( 6 « l) 


so timt 


For this fanction of y > the effective dispersion coefficient is 
calculated as 

]y{- ^ G*(a,g,5) (4 34) 


where 


G*(a,0,6) = = ^ (4 35) 

( a sinh a + 6 cosh a ) 

A* =-U-l{a^ - (a^-f3)(l5a^ + 5a^B + 136)) sinh a 
45 a 

+ (45“^ + 150^8+456 + a^8 ) acosha] (4 36) 

a nd 

- 6 (“V ^ ° 

60 ^ 6a^ a' a a 

+ e(~4 + ^ + i-)] (4,37) 

36a 3a 3 a a 

When fi = 0, Gf*(a,e,S) reduces to 


( 5 *(a j 6) rs -1- (a* suah o + B* cosh a ) 

I cn ‘ * 


(4 38) 
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^riiere 

a® ^ (a^ +5) (i 5 a^ + 5a^8 + 15 6) 

^ 4-3]^ (P( a sxnb a + 6 codi a) 

0 2 4 

45^)t + "j5 ^ B 45 0 + ct B 


B* 


^ ' 

^ 45 V a {a sioh a + 0 codi ct ) 


(4 39) 


(4 40) 


vih.xdh can be seen to agree wxth Gupta et al (l972) 


As 0 0, we get 


/ n 1 r /1 1 1 coth a 

It G*(«,0,fi) ='^ t (45 - 2 ' 4 ^ 

3 -*-0 “ 5a a or 


^56 ' 2„2 - „4 


4 cosech g 


2 colii a , 4 coth g 




01 


(4 4l) 


winch agrees with eg,uatien (4 24 ) 


2b see Ihe effects of 8 and 5 on Ihe effective daspfflpsioa 
coefficient, -fee function G*(a,8,6) given by equation (4 35) is 
calculated for 0= 1 0 and tabulated in table (4 2) 

It can be seen from table (4,2) Ihat, for a given a , G«<a,Bt5) 
decreases as 0 increases ifdiereas it increases with increasing 5 
when 3 is held fixed 
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Sable 4 2 

(a= 1 0) 

8 

5 

G*Ca,B,6) 

2 0 

0 00 

0 16969080 X 10"^ 


0 01 

0 17172264 X lO”^ 


0 1 

-2 

0 19000921 ^ 10 


0 2 

-2 

0 21032761 X 10 


0 3 

-2 

0 23064602 x 10 


0 4 

0 25096442 x lO"^ 

4 0 

0 0 

0 16605727 x 10“^ 


0 01 

0 16804359 X 10"^ 


0 1 

0 18592C50 X lO"^ 


0 2 

0 20578372 x 10"^ 


0 3 

0 22564695 lO"^ 


0 4 

-2 

0 24551018 X 10 

6 0 

0 00 

0 16457162 X 


0 01 

0 18053926 X ld"^ 


0 1 

0 18424333 lO"^ 


0 2 

-2 

0 20592444 10 


0 5 

-2 

0 22360085 X 10 


0 4 

0 24327726 x lO”^ 



1CX) 


8,0 


10 00 


00 

0 16376549 

-2 

10 

01 

0 16572301 x 

lO'^ 



-2 

1 

0 18334072 x 

10 

2 

0 20291595 X 

_2 

10 

3 

0 22249119 X 

10“^ 

4 

0.24206642 x 

io"^ 

00 

0 16325898 x 

10"^ 

01 

0 16521009 =x 

,o-= 

1 

0.18277007 X 

10"^ 

2 

0 20228116 X 

,0-" 

5 

0 22179225 x 

id"^ 



-2 

4 

0 2413C934 X 

10 


HSjos it may be concluded that the effective dieperslon 
ooefflciettt decreaaes due to hosnogerteous and heterogeneous 
reactions and xt increases as the riscosllgr of the flixid decreases 
tovrards the 'ssall 

4 4 ErTECfS OP HOhOGEHBOUS RZAGI102J py DISPIS'iION A HilTlD 
ELO^dlTO lU A CIRCULAR IDBE 

In the last two sectior® we have studied the effects of 
hcMogeneous and heterogeneous reactions on the dispersion of a 
solute m a fluid flowing between two parallel plates Ctorrespondxi^ 
results will be obtained in this section when 1be fluid is flowing 
in a circular tube. 


I 
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Oonsider the laminar flow of a viscous liquid of -vaiying 
Tiscosi-ty under constant pressure gradient in a circular pipe 
^ose physical configuration is sitown in Eig (4 2) is befctre 
im consider the visoosily variation as a known function. 


p =y(r) (O £ r <_ R ) 


(4 43) 


Ihe fully developed velocil^ profile cem he obtained in this 
case as 

- 2 

r ^ 

Hie average velocity, #ileh is defined as 
E 

W = ^ J r?dr (4 45 ) 

E o 


is given by 



¥e assiaae tiat the solute diffuses wiule undeigolng 
reaction in the liquid under isothermal conditiom 
for concentration c is given by 


(4 46) 

irreversible 
Hie equation 


2 L. 

r 3r 


(r||) -Ko 


9c 


+ JC 


(4.47) 


Relative to a plane moving with, the mean ^eed of the flow, the 
fluid velocity is given 




f(n)l 


X 


(4 48) 




lIPMi 




vmm 




X ,v;;' '/ ■■! 
' ■/ ■';: : 
muiBiniii 


IIP 

sill 




XX0. ‘XXXp i. 


in 
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wbere and fC^l) are given by equatxons (4 7) and (4 s) 
Ifondxmensoonalxsxng the variables by 



n 


r 

R 


and following the dlaylor's ^preach equation (4 47) can be 
written as, 


1 £- (n ^ ) _ a^c = H [G - f (n ) ] 
n 3n ° 

where 

“ D 

Ibe boundary oonditions are 


(4 49) 

(4 50-a) 
(4 50wb) 


=x 0 at n=: 0 and n «= 1 (4 5 1 ) 

3n 

SolTring equation (4 49 ) with 1die abore boundajy condition, we 

get 

c = [fro] Iq («n ) - — ~ - H p(ti ) (4 52) 

i.^v / 

^diere 

S =s / t K.(a) + I^(a) K^(ot)J f(t)at (4 53) 

o ^ 

T) 

P(n) = / t [I^(an) Kjert) - I^Cat) KQ(an)l f(t)dt 

° (4 54) 

I . K are aero order modified Bessel functiojis and I. , K- aio 
0*0 ' 

of order one of first and second kind respectively Ihe 

■y 



1C4 


volumetric rate at wiixcb. liie eolute is transported across a section of 
Hie tube xs 

E 

Q = 2ir J r c W dr (4 55) 

o 

Using equations (4 48) and (4 52) we get 


,8 


« = - ir (if) 4») 


(4 56) 


wtiere 


1 SI (an) 

x(a) = / nCF(n) - - 3 [% - f(n)] 


dn 


(4 51 ) 


o 1 

Oomparing wiHi Fick's law of diffusion tie effective dispersion 
coefficient is obtained as 

DK- = (■^)^ “ X( a) 

2 '■dz'' D ^ ^ 

where x(a) is given Igr equation (4 57) 


(4 58) 


If Hie fluid has a constant viscositXj expressions for concentration 
and equivalent dispeimon coefficient can be obtained, after scaae 
manipulation with Bessel functions, frcaa liae equations (4 52) ard 
(4 58), as 


% o’ 1,(0) ° ai^ / ^ 0^ 


(4 59 ) 


and 


52 L(a) ^ 

Y^en o is Email, Bx- can be written as 
48 15 ^ 


(4 60 ) 


(4 61 ) 
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We observe ifoai a= 0, the Taylor’s result xs obtained We rsotice 

* 

also as ct increases D decreases 

It may be noted that the e3q)ressxon for D* given by 
equation (4 58) xs for any known function of viscosity variation 
(fenoted by yCn) However, to study specifically the effects of 
viscosity variation on dispersion the following relation is assumed 
as before , 

W (r) = (6<<1) (4 62) 

or 

j^ = ^(l+5r) (4 63) 

50ae effective dispersion coefficient is calculated from equation 
(4 58) using the equation (4 63 ) for viscosity variation. 

!lb see the effects of a and d on the effective dispersion 
coefficient, polynomial expansions of the Bessel functions have 
been ^^sed to evaluate the integi^ls in x(o,d) py neglecting 
and higher powers The values of ltx:(a,d) lave then been calculated 
and tahxxlated in table (4 3) 



106 


3^1e 4 3 


6 


0 00 

0 65060771 X ICT^ 

0 01 

0 65924755 x 10"^ 

0 1 

0 74129019 10”^ 

0 2 

0 84149353 x 10"^ 

0 3 

0 95121762 X lO"'^ 

0 4 

0 10704625 lO”^ 

0 00 

0 64930735 x icT^ 

0 01 

0 65793009 X kT^ 

0 1 

0 73981530 10"^ 

0 2 

0,83983510 X 10”*^ 

0 3 

0 94936667 ICT^ 

0 4 

0 10684100 X kT^ 

0.00 

0.64714460 x io"^ 

0 01 

0 65573^8 X id"^ 

0 1 

0 75736230 X 1<r^ 

0.2 

0 83707686 x io"^ 

0 5 

0 94628819 x fO"^ 

0 4 

0 1C«49963 X ICT^ 


0.3 
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a 




0 4 


o»3 


0 00 

0 

54412641 

X 

10-^ 

0 01 

0 

65268124 

X 


0 1 

0 

73393912 

X 

ICT^ 

0 2 

0 

83322773 

X 

tO-5 

0 3 

0 

94199227 

X 

10"^ 

0 4 

0 

10602327 

X 

kt^ 

0 00 

0 

64026262 

X 

10-5 

0 01 

0 

64876683 

X 

10-5 

0 1 

0 

72955689 

X 

10-5 

0 2 

0 

82830031 

X 

icr5 

0 3 

0 

93649294 

X 

IQ-^ 

0 4 

0 

1C641547 

X 

icT^ 


It caa be seen fixaa table (4 3) that as a increases X(a, 6) 
aecreases for a fa:^d 5 whereas it increases as s increases when a 
18 held fixsd, 

4.5 DISPIP^IO^ 0? A 33]3JfE IS A mril) gm^r^G A TIPE ffHH (X71BPE'^ 
H0''^0SEI7’X)US AND HFIIROGE^rOUS FEACTIOi 

In "tills section we study 1iie dispersion process by t^ing 
into account the effects of chmical reaction in lii© fluid and at "Sie 
wall of the tube 

Ihe diffusion equation in this case r®iains the same as 


before, i e , 




1 (» 


(n|^) - c = H [G - f(n)] 

n 3n ° 

(4 64 ) 

She boundary conditions beccaae 


■2E = 0 n « 0 

9n 

(4 65) 

■^+gc = 0at n=1 

9n 

(4 66) 


where B = sR aid s is the reaction rate constant 

Solving equation (4 64) and using the boundaiy conditions, (4 65) an«i 
(4 66), we get 


H G 

c = A IgCoh ) ^ - H P(n) 

a 


(4 67 ) 


where H and ]P(n) are given by eqmtions (4 50-b) aM (4 54) 

respeotavely and A is given by 
G 


h[B — ^ + 6 P(l) + aS] 


A = 






1(1) = / [IqC^) E^C^t) - I>t) K^(a)] t f(t) dt 


(4 68) 

(4 69) 


and S is given by equation (4.53 ) 

tphpt flux of the solute across a section of the tube is given Igr 


8 


n _■:££- Y(a.B) 

^ ■■ ZHi ^dz^ ^95'^ ^ t J 


(4 70) 


^er0 


Y(a,6) ^\ict) - Y^(a,8) 


(4 71) 
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0-^ + BF(i) + aS ^ 

Yi(a,B)=[ + g i^(oj I / n Eg-^ - fCn)] ^aT,)d n (4 72) 


a-id 


1 

y (a) = / p(ri) [g _ f (n )J n dn 
^ o ° 


(4 75) 


Gcmparing with Pick's law of diffusion, the effective dispea^sion. 
coefficient is given by 


Die- = S— (j^y 
21) Mz'' 


Y(a ,S ) 


(4 74) 


In order to study specifically the effect of viscosity 
variation on dispersion in the presence of homogeneous and heterogeneous 
reactions, re assunie that the viscosity variation is given by equation 
(4 63 ) Polynomial eicpansions of the Bessel functions have been used 
to evaluate the integrals m Y(a,e,6) by newest ing a and hi^er 
powers Dae values of l(a,8,6) have then been calctilated and tabulated 
in table (4 4) 
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(fefele 4 4 

(a = 0 2) 

6 

f(a,e,6)= 

0 00 

0 64717996 x icT^ 

0 01 

0 65377342 ICT^ 

0 1 

0 73738613 ICT^ 

0 2 

0 83708411 1<r^ 

0 3 

0 94627387 ICT^ 

0 4 

0 10649554 ICT^ 

0 00 

-3 

0 64716943 x 10 

0 01 

0 65576275 x kT^ 

0 1 

0 73737410 X ICT^ 

0 2 

0 83707049 X 1(f ^ 

0 3 

0 94625^5 X kT^ 

0 4 

0 10649382 X iCf^ 

0 CO 

0 64716586 10*"^ 

0 01 

0 65575915 x ICT^ 

0 1 

0 73737005 x icT^ 

0 2 

0 85706590 x 10"^ 

0 3 

0 94625341 X 

0 4 

0 10649325 X ICT^ 



111 


3 


8 0 


10 0 


6 

Y(a,e,6) 


0 00 

0 6A716411 

X icT 

0 01 

0 65573736 

X icT 

0 1 

0 737368C5 

X 10“ 

0 2 

0 837C6362 

X icT 

0 3 

0 94625084 

X io" 

0 4 

0 10649296 

X iCT 

0 00 

0 64716306 

X icT 

0 01 

0 65575627 

X iCT 

0 1 

0 73736683 

X icT 

0 2 

0 83706222 

X icT 

0 3 

0 94624929 

X id" 

0 4 

0 10649279 

xicf 


It can be seen froaa table (4 4) ttiat Uie effectxTe 
dispersion coesffxcxent decreases as P increases whereas it incrmaes 
with xncreasit^ 6 for a frxed 0 Similar result has been obtained 
in Ihe case csf paiallel plates. 
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6 (Dlicnjwl I tS 

Eic effect of viscosity v'^riation or luylor diffusion has 
Decn inyestigT,ted in tiie cases of flow between parallel plates ard 
circular pipes by taking onto consideration tbe bcoogeneous and 
leterogenoous reactions of the solute with the fluid. It is found 
in the above cases that the decrvas>- of viscosily towards tne vsall 
increases liie effective dispersion coefficient 

It IS also observed that for a given viaiosily distribution, 
the effective dispeieion coeff ic lent decreases with 1iie increasing 
reaction r^te constants withm the fluid and at Ihe wall 





Goiicertration of tne solute in the fluid 
oleeular diffusion coefficient 
Taylor’s equivalent dispersion coefficient 
Half of the width oi the channel 
Reaction rate constant in the fluid 

HLux of the solute across any section of the diam^ or tube 

Radius of the circular tube 

Telocity of the fluid between parallel plates 

Average velocity of the fluid in tte parallel plate case 

Relative velocity of the fluid in parallel plate case with 

respect to a plan© moving with the mean speed of the flow 

Telocity of the fluid in circular tube 

Average velocily of the fluid in tl» circular tube 

Telocity of the fluid in circtilar case relative to a plane 

moving with the mean speed of the flow 

Ibn-difflensional reaction rate constant in the fluid 

Kin-dimensiODa]Used reaction late constant at the wall 

Tisoc»itv of the fluid 

Tisoosity of the fl^^ld at the central line 
ISfon-dhaensionalised asiial distance 


Fon-duaenslonalised transverse distance 
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DlSPiPSIOh IH CASSDlf "^ID 


5 1 IlTOOHJCflOI 

In Chapter If, we have studied the dispersion iirrough a jiiysio- 
logical lewtonian fluid with variable viscosity, flowirg throu^ a 
rectangular chanml and a circular pipe In general, biological fluids 
such as blood, do rsDt behave as Newtonian fluids [ Colcelet et al 
( 1963 ), Meirill et al (l963), Cham atd Zurland (1962 , 1965)1 and 
can be represented by a non-Newtoman Casson model fluid at low 
shear rates* (See CSiapter l) 

As pointed out in Chapter I? the dispersion thiough non-lewtonian 
fluids flowihg throu^ ducts have have been studied by different 
investigators by considering different fluid models, power-law fluid, 
item and Hwang (1965), Bin^aa and Ellis models, Fan and Wang (l966), 
Beiner philippctff model, GShodhal (l97l), E^ing model, Siah and Cox 
(1974) 

Eeepii^ this m view, in this (Siapter we study toe dispeirsion 
of soluble matter in a Casson model fluid, followirg toe approach of 
Saylor (l953) 

5 2 BlESPEESION IH CAS80H MOISL FIDID EliOWINg BEBifEM WO PLAfES 


Consider the laminar flow of casson modal fluid under a constemt 
pressure gradient between two parallel plates distant 2b apart, toe 


(17 


physical confi^ratxon aM the coordinate systeaa are ^own in 
fig. (5«l) 


Bie filiation of Biot ion for one dimensional steady Isainar flow 
of Casson model fluid is 


-Sel „ ii = 0 

dx dy 


(5 1 ) 


IThe atress-R^aia law as given bj 



=y(- 


dV,,l/2 







(5 2 ) 


Where t is the yield stress of the fluid. Integrating equation 
0 

(5 1 ) and ttBxrg the boundary condition 


T =s 0 at y *= 0 


we get 


T a? 


(- ■^) y 


(5 3 ) 


If y^ is IdK hel#t of Idle core then we have 




(5 4 ) 


Eliminating t from equations (5 2) and (5 3) and integrating with 
no slip condition at the wall, we get 

V . 1 [£ (iV) + T„(h.y) - I ^ (5 5 ) 


2 ‘2 




'a ' 3 * ^ 

" ?‘ . ”' U' tif ‘' iPff I 'If ’ : 

'vwC ^y‘,iKV 




mWM 




SMh^b 




!V‘V; .i""'//'?^". ; 




';-yv-;/y2' 




llil 


*iii 
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in the legien (y„ £ y S b) ^ere I' = - ^ ®ie velocity 


co 3 ~e is 03 -'veh by 




in the region (o <,y l.y^) 


Substituting 


n = y/h , tIq = yyb 

in equations (5 5) and (5 6) and dividing -m obtain 
^ (i-n^) + 2n ^ (i-n) - ^ (i-i?'^^) 




C^io i ^ 5. 


whOTe 


8 1/2 1 2 

1 ->■ 2t\j - 5- - 3 % 


and 


y 

^ 1 

0 


(0 5b < n ^) 


Ihe average velocity f , #xich is defined by 


IS obtained as 


y « m Y„ 


1 , , t 

Sh 


-h 


y dy 


If* ere 


_ 1_ / 2 
» « V 3 + be 


24 1/2 

15 % 


3 X 

- 15 bo ) 


of the 

(5 6) 

(5 7) 

(5 8) 

(5 9) 

(5.10) 

(5 11) 

(5 12) 


(5.13) 
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’Bie eq.'uatioii @0Te3?iiing "the coocen.'fera'tioii o of the solute is given 
by 


3c 

■^t 


+ V 


3c 

3x 


2 

ay 


(5 14) 


Relative to a plane laoving with the mean speed of the flow, the 
diffusion equation becomes 


-^ + (V _ V ) 4^ 

at: -'ax 


Hbndnaensionalising by 


D 


ifs. 

ay^ 


(5 15) 


x-Vt 


t T 


, e = ^ , t =r 
t V 


ard follx)wing ^proach the dlffuaion equation becomes 


2 2 

l-£ a iL (M.) V 

3^2 m X 


(5 16) 


where ^ ^ is assumed constant, 

yrcw. equataons (5 8), (5 10) and (5 12), can be obtained as 


follows 


Y = V ^ ^ 

X 


(n^-h^) + 2h^j(TWTio) - I t^^^) 


m 


So 


(5,17) 


in region (n^ i h 1. 1 ) 


= ^ (i- 1) (5 18) 

in region (O h ) 

¥4iei« g^ and m are given by equations (5 9) and (5 13 ) 



121 


Bie boundary cormtioas are 


(x) 


ac 

an 


= 0 


at n = 0 


(5 19a) 


(xi) = 0 at n = 1 

an 


(ili) = c at n= 0 


(5 19b) 
(5 I9e) 

(it) c xs continuous at n = Hq ^5 l9d) 

SolTxng eq,uation (5 16) usxi^ eqtiatxons (5 1?) through (5 I9d), 


we get the concentratxon c of the solute as follows 

“ = vs 

(ilo 1 1 1 1 ) 

[(-- 1)%-1 + c 


(5 20 ) 


m 95 


m ■' 2 ■ '0 

(o < n iho) 


(5 21 ) 


where 


fin; = - “y~ - 105 % - 6 ^0 ^ 


1 3 1 4 

+ n-^ n^ -sr nZ 

15 o 84 o 


(5 22 ) 


TOluaetric flow rate at ^xch lire solute is transported across 
a section, of the channel df unit breadth is 


% 1 
2j oV dTi+2 / cV^dn 

o n_ 


(5 23) 


Substituting appropriate expreseions for c and aod Integratxng, 
we obtain 
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. 2h^ I 

^ = ^85^2 I 6 




where 


En 

o 

84- 

153 


7E 

En 

_£ + 

344 E 1 

' 60 " 

5 ^ 

945 

+ J_ 

n 

+ -2. 

332 

84 

12 

" 4095 

I 0 

0^2 

15 

944 
“ 3465 


l/2 ^ 3 

I E n^ 

_o o 

30 

.3 


60 


,3/2 


o 

84 


+ I22n -i- . — 

945 o 315 o 315 o 6 


l2n3 


29B 1363 tlo 


270 


77220 


and 


1 l6 1/2 1 3 

+15% 


(5 24) 


(5 25) 
(5 26) 


Goapam® with Hok’s law of diffosioia, the effhctiTe dispersion 


ooefflolenfe is gxyen by 



2 -2 

h r 


D m 




(5 27) 


where 


J<\) 





When n = 0, the effeotiYe 
O ' 


dxspersaon coefficient is given 



(5 28) 
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iwbxcii agrees with the result of Wooding (‘i960) If << 1, then 


h and hi^er powers of n are saaLl when ccmpared to h 


,1/2 


Mth this ^prozxsHtion, P xs approx inated as 
.2ri 


D' 


.* _2_. 128 _l/2 


n " “ ^ 
195 o ^ 


(5 29) 


I> 105 

It may be oo Deluded fresa equation (5 29) lisat the ToruHewtonian 
character of the fluid decreases the effective dxspezmon 
ooeffi'^ient 


lo stu^ the effect of non-Hewtonxan character on effective 
aosperstor. ooeffioaeirt, X ^ oal,mlated and talmlated In 

M 

table (5 1 )• It can be seen from t^le (5 1 ) that as imreases 
the effective dispersion coefficient decreases 


5 5 hlSPlT^SIOir Itr CjiSSOH ^OICL -ILI!) Sih p^r^iCE Or 
HO oullESUb PEA^TIQj. 

In ihis case we assame "Siat the solut f while diffusing, 
undergoes an irreversible chesaicaL reaction in the fluia, Ihe 
diffusion eq\xation then becemes 

ff+ (5 30) 

ay 

Poliowing Ifaylor’s appjroaoh, as in the last section, the diff- 
erential equation governing the distiabution of the soluble matt- 
er relative to a plane moving with the mean s^ieed of the fluid 


is obtained as 
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3 c 
3 n 


2 “ KL ^ 35 X 


(5 3l) 


v^iere 


^2 ^ Kh 


B 


(5 32 


and is given by equations (5 1?) and (5 18) m +hs two 
regxor® 

Itie boundaiy conditions are 




at n = 0 (5 33a) 

n = 1 (5 33b) 

Q (5 33c) 

Solving liie equation (5 3l) witti the boundary condition (5 33a) to (5 33c) 

we obtain the cwancentration of Ifae solute as follow 


(li) 1^=0 

(xii) c am |£ are cortlmoue at n = 


c »s H [ig cosh on + Bg sirih an + P(h ) 1 


< n < i) 

o — 


(5 34) 


lb ere. 


A„ = t ,. -- (m - p oosb an„ oo^ “ + ®xnh an . cosh a ) 

^ sub a 0 o 0 o ° 

(5 35) 
(5 36) 


Bg = p^ oosh <*n^ - sdjb Wq 

H = ^ V ( 1 ^) 

HL '■ 35 ^ 


(5 37) 


- 8 n 


, 1/2 


aM =s 


0 3 n n 


1 ,/^ 2(1 + ru) 

/ xr co^ (ou ax) dx - •' —" " g (5 38) 
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8n‘ 


3a m m a gp 

4n 


o 


3 

mag 


(5 39) 
(5 40) 


f(^) 


a»Dd 


1 2 / ^ Jf ^ ^ T 

= — [ 3 “ (n + n ) + 6 _ 4 anQ ] 

3h“So 

1^2 

" “2 ^ ■' 

a 5 iaag^ 


r = H cosh an- % - 1 ) (0 < n 1 n^) 

‘ a 


(5 4l) 

(5 42 ) 


where 


A^= cosecha- cosechan^ 

-fp cosh on - H sinh an„) (coth a - coth an„) (5 43) 

o 00 o ^ 

and H ,P .H are sane as in equations (5 38) to (5 40 ) 
o' o 

Wien ri = 0, TO get the conceiatration c as 

o' 

” = -2 (|f) (|f) * -^|f)(f - % -"^3 

mr V sirii a 2 v£!L o a (5 4.4 ) 

which ag3?ees with the result of Gupta et al. (l972) 

Tbe vohMetrdo flow mte Q at whidi the ecaute is transported across 
a section of the ctonnel of tmit breadth is given by the equation (5 23) 

aibetituting the Expropriate expressions for o and from eq\iations 
(5 34 ) and (5 42 ), (5 17 ) and (S is) we obtain the flux Q as follows 


p(a,n ) 

^ Ml o 


#iere 


(l-m)%Q A (l-m) sinh oHo j- 

p^a^n^) _ 2 2 “ ' 


m ct 


a 


2 2 2 
m a g^ 


{5 44 ) 


(5 46 ) 


j is gLven by the following equation, bl £ind 
in equations (5»13) and (5 43)* 


are same as 
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J =(8illh «){— . 




a 


J r.2 

v2 


,, , .,'>l'2 =''+"0^ t’+'t) =2 ^"2, 

+(oosh a){— - + r _ — } 


8 l/2 

+ 5 V 


a^ “a? 

/ (£, coHh cw =2 sliil a. 


fti^^ /’ (♦,-ti+n^)^ + I T(n)aii 


,<> e 4n e 2e 4n 

.(311* c.nj(— -—^- — + -^-) 

a a 

. (cosh an^ + — 2 — ~T~ ^ 

a „2 a „5 


8 nVa r 2 2 g 

3 o h W 5 


4n 


o . 2 ^2 


+ V n:' ■' l^ + # + ~ 


4 e 


- n^ (2 -If ^■+€ 5 ) + tIq (3<I>^- 2 + ^ e^) 


aod 


4 1849 5 

*■ ’'o " 189 ’’o 


e tta iSq 

2 ^ 

e = ma g B. 

2 o 2 

e = ~ - — - (iJa) g^ 

“ >i«2 

(l-m) g^ + — 

n sihb a(n-x} 

®(n) = / 

n. ® 




dx 


(5.47) 

(5 48a) 

(5 48b) 
(5 48c) 
(5 48d) 

(5 48e) 
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Comparing witn Pidc’s law of diffusion the e ul-yalent oispersion 
coefficient is given by 

(5 49) 

where is giv^n by equation (5 46) 

When = 0, frraa equation (5 49),D* can be obtained as 
follows 


* __ fh^¥^ 1 c _1__ 1 coth g 1 1 

^n=0 “ D 2 '■ 45 " ,_2 5 " 4 ^ 

a 3® a. a 

which agrees with the remit of Gupta et al (l972) 

Iben is SEsall such that and higher powers are 

1/2 

neglected in comparison to , the eqiiivalmt dispersion 

coefficient raiuces to 


(5 50) 



h^ 

D 



(5 51) 


where 


V (a n \ rri- -i~ + ^ ° L.\ 

2 2 2 ^H5 ’■_2'*' 3 “45/ 

m s g^ 3a a 


0 ^ 

3 5 a2 


2S coth g 23? coseoh a 
a a 




</)] 


(5 52) 


and 


/ t5/2 ooeh (o -at) dt 


(5 53) 
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T 


= / coati (at) at (5.54) 

*0 

= aia. (»t.as) ^3/2 (555J 


O 


As a •* 0, it cm be i*o?m that 


Lt D 


b^ 2 /, 128„1/2 n 

(1 » “ ) 


1 D 105 


a -► 0 

wbich xs sane as eqmtion (5 29) 


195 ^0 


(5 56) 


To study the effects of homogeneouB reaction on dispersion, 
equation (5 5l) is calculated and M>ulated in table (5 2) It 
can be seen frcoi table (5 2) Ifcat as a increases, the equivalent 
dispersion coefficient decreases md "Bus decrease is ertonced 
furlijser by increasing 

5.4 Bisaernsrog ii gassos hpib Hoim siRcoffl i ciamAR 

(KIBE 


In Idle last t«?o sections vre have studied the dispersion of 
a solute in a casson model fluid fLoidng between two paiellel 
plates In ^e foUowing we get correi^nding results when the 
fluid is flowing through a ciarcailar pipe She fully dereloped 
velocity uTOfile in Ihis case is given by, (see Eig 5 2)) 

7 = -L [| (E^-r^) + T^(E-r) « r^^^) 1 (5 57) 


3 2 


in -the region (r^ ^ r <.E) where ^ ~ ^0 

radius of the core» Whitmore (l968) 




wmMm 




S , C "4f‘^>V Ai^ 


'W'if ?,if;‘'i 






150 


Bie relocxly of thB core xn the region (o 5. r £. r^) is 
gxTen hy, 


~ [f Tq) + - 5^2^® ” *0 


,^3/2 _ 3/2 


)] 


(5 58 ) 


In this case, the yield stress 

oore are reOsted by the relation. 



aobstituting. 


t ard "ttie r^idius of the 

o 

Bird et al (i960) 


(5 59 ) 


n = r/E , rig = and dXTXdir^ Ihe equation 

equation (5 58), we get 



in the region (n^ ^ 


(5 57 ) by 


(5 60 ) 


I 1 (5 61 ) 

f 

0 

in liie region (O in in^), where is given by equation (5 9) 

®ia average velocity, which is defined by, 


1 

1 = 2/ nw dn 

o 

is obtained as 

W - n Mr„ 
c 


(5 €2) 


(5 63) 


where 


n 


I? + ” 7 “ - 7 





(5 64) 
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rollowiig Ifeylor's approach (l953)i "Kie govemxng etiuation for 
comentration of the solute relative to a plane iMving mth tlte 
mean speed of the fluid is given by 


3^c , 1 ac 

2 n ^ 
9n 


35 


w 


(5 65) 


Tstoere W-W) is the velocity of the fluid relative to toe 

moving plane and is given by equations (5 60), (5 6l) and (5 65) 
as follows 


= l) - 


n 




(n^+ 


8n 


1/2 


2Ti^n- 








(5 66 ) 


and 

W = ¥ (-J- -i) (O 1 n £n^) 

X U ' > ^ Qf 


(5 67) 


file boundary conditions are 


(i) ~ fs 0 at n= 0 (5 68a) 

(li) 1^ s= 0 at Ti!= 1 (5-6^)) 

(iii) e = Cq at n » 0 (5 68c) 

(iv) 0 is oontiiiLsous at pstJq (5 68d) 


Solvii^ the equation (5*65) with the bmindaiy conditions (5 6^) 
through (5 68d) we obtain the concentration of toe solute as 
follows 


(5 69) 
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wQere 


g(n) =ji + !^.32_„l/2 „7/2 

16 9 147 ^0 


ii5n 


^2 ^2 ^4 

12 84 + 7® 6 


and 


tj „ ^ i® 
H. 35 


c >. 5- ^ If [(-i^S) —I + c (o<n<n) 


(5 70) 
(5 71) 

(5 72) 


13ie Toluiae'fecic flow rate with ?diicn the solute is trans— 
poirted across a section of the tube is 


2 ^ 

2 irH J n c dn 
o 


Substituting the expreesiom Ibr c and 
(5 69), (5 72), (5 66) aM (5 6?) we get 

^ ^ ^16 "o ^ 


(5 73) 


froB equations 


(5 74) 


irtiere 


® ^1^2 ** 2^^ " ^Vo ” 3 ^ 


,4 


+ (44'^ - 4^) - -f’) 


26ti 




4nf 54, 
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(1 

1176 o 7 


7 8 
-o ) + (_1 2 ) 

7 '■128 128^ 


5 0 ° ^7 7 % ^ 


8 l /2 \ ^ J 1 / 2 V 


t5€0 3/2 r2 2 13/2% .115_ 1/2 15-^2% 

- I - SIO * ”n 1 ^ 1 ' Hn vl “Ho ' 


1325 'o M3 13 
^ o \ 

-—g^c- 1 - log nQ+ r,^) 

5 n^ 2 

. % ( 1 '’o , „ ^ % % 

1 --^C - 4 - -"T ^ 


228 $ ° 


2 n ' 


.9/2 


63 


\— 2^ 


4n 


2 5 / 2 .^„ . - 

“ 5 ’’o + -^ 


3/2 


.4 


■*■ *51 C- ^ log n„ + *— ) 


84 


(5.75) 


where 


84 

■ 3 w 

8 


log n,- 


115 4 


b 7056 


8 1/2 
2 T’’o 


+ 'r“ - -rmr + *755 n. 


1 4 

168 *^0 


(5 76) 
(5 77) 


CkaHj^riDg with Slefc’s law of fiiffiaaxoQ, the effective 
dispersion coeffioieKt is obtaired as 

= ^%*C*A-) (5 78) 

8n ” 


Q g . 


When « 0, we get the sff ectire dispersion coefficient 


as 



154 


D 




f± 
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winxsa agrees with the result of !I!^lor (1953) 


(5 79) 


ibr « 1 aod neglectirig rip il^s higher powers ia 

A ^ 

ccMpariaon to , D can be approxonated as 


W 


r t . 5g_„l/2i 

1 '• 48 " 5465 o ^ 


(5 80 ) 


Hence, it may be comluded in Ihis case also that iiie ron-hewtonian 
character of the flow decreases the oCfwtlta dispersion coefficient 


So study tne effects of Eoo-Hewtonian bdiawiour on dispersion, 
equation is calculated and tabulated in table (5. 5)* I't 

can also be seen from the table (5 3) ttot as irKireases, the 

dispersion coefficient decreases 


5 5 igsggtqiOF IF C4S*^H POtElt fPnP) 'PTOVire IB A QISGOItffl OTBB 
I 7 3HE g-, - ■ OE Q-t (HB-ilCj-vL FlbiCTIOIT 

In this section we study the dispersion through a oasson fluid 
by consideriEg tiis eff^ts of an Irrerewible diaaioal reaction 
in 1b® flttldU 

aie differential equation gorerning the oomentietion xs 
aiodified as follows* 

^ + (0<i,<n„) C5 81 ) 

3n 

+ i is, „ o^Q * H r (i - i) - 5. ^ 1 0 (5 82) 

2 n an ^ ag -i o ~ 

an 
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where 



&(n) = n^ 

8 2 3/212 
+ 2r^n-jnQ n - 3 no 

(5 85) 

She boundary conditions are 


(i) 

Is. n 

1! 

0 

(5 84a) 

(ii) 

.^=0 

3n 

at n = 1 

(5 8lb) 

(iii) 

3c 

c and - 5 ^ 

are oontixiaoas 



at n 


(5 84c) 


SolTirg the equation (5 81 ) and {5 B2) and usir^ Ihe boundary 


conditiooB (5 84a) to (5 84c), we obtain the concentration ae 


follows 


o 


0 


H S^I^(an) 

hg^X^C") 


B (vn) 

s— (0 < h in ) 

nor 


S( Vn) 
no 


-S- f (i*) 

% o' 


(h in < i) 


(5 85) 

(5 86 ) 


S = f (at) + I (at) 1 (a)3 tG(t)at (5 8?) 

1 1 O U I 

F (n) » / [l^(un) K^(at)-lQ(ot) K^(an)ltG(t) dt (5 88) 

° 

and H and n are g^r^a by equations (5 7l) and (5 64) 


respcKstiTely 
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She iroltmetric rate at viiidi tfle solute is transported acaxias 
ai^r section of the tube is given by equation {5 75) axbstituting 
appropriate eaqjressions for e and W from. eq\iat±ons(5 (5 

(5 66) and C5«67) we obtain Q as 


=ss!£ (1£) T (c.,n ) 


{5 89) 


Titiere 


E(t) io(«) it 

(t-n) 


+ (l-n) go / 


/ G(t) dt 1 


(5 90) 


Goaparing wildi lidfe's law of diffusion, tl» offeetiTa dispersion 
coefficient is obtained as 

^ Y (a,n ) (5 91 ) 

B o ' o 

whftw n t:; 0, we obtain Idle effective diversion ooeffident 

-A * 


frc» ©qmtion (5 9l) as 


B 




eV j32 IjC-) 


] 


(5 92) 


“ 11=0 D ‘^5 5^2 

is sane as obtained in Chapter IT, equatioa (4 6C) 

mien o-*- 0, and Bq is aaall so that and hl^r powers 
can be ne^eoted in coaparison to , we obtain Ifae effective 


f 
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diepersxon coeffixiient as 


a=0 " D ^8 “ 5465 ^ 


wiilch. is the stsae as in equation (5 


(5 93) 


©3 stud/ the effects of hcmogeneous inaction on dispersion, 
the equation (5 90 ) is calculated and tabulated in table (5 4) 

It can be seen frcsa table (5 4) that as a increases, the effeetiTe 
dispersion coefficient decreases and this decrease is enhanced 
furiher by the increasing 


5 6 OOHCKFSIOBS 

In this chapter, dispersion of a soluble aatter throu^ a casson 
atodel fluid flowing in a rectangular channel and circular pipe has 
been studied 

It is fotiM that the effectiTe dispersion coefficient in ti^ 
oasson fluid is less liian that of lewtonian fluid, ^en = 0» the 
xeailts are found to agjjee mth those of Gupta et. al {l^2), 

W>odihg (i960) and S^lor (l^3) It is aim obserroi that in the 
presence of hcMOi^necRis reaction, Ihe equivalent diapersion coeffioient 
decreases when the reaction rate constant increases This decrease is 
enhanced further by the increasing non-Bewtonian character of the 
fluid for a given sean v^ocity 



Table 5 1 


Ho 

■wr o' 


0 00 

0 19047619 X 

10"^ 

0 02 

0 17070296 X 

10"^ 

0 04 

0 16116019 

10"^ 

0 06 

0 15324750 X 

10"^ 

0 C» 

0 14616516 X 

10""* 

0 1 

0 13960928 

10'^ 

0 2 

0 11110413 X 

10*“^ 

0 5 

0 86673169 

10"^ 

0 4 

0 65135510 X 

-2 

10 

0 5 

0 46335787 x 

-2 

10 

0 6 

0.30405961 * 

10 



Table 5 2 


% 

a 

%Ca,n^) 

0 02 

0 1 

0 11239^6 X io"^ 


0 2 

0 17201045 X 10"^ 


0 5 

0 17103077 X 10"^ 


0 4 

0 16965921 X 10"^ 

0 04 

0 1 

-1 

0 I6527I88 X 10 


0 2 

0 16467729 * 10“^ 


0 3 

0 16368629 X 10”^ 


04 

0 16229890 X 10“^ 

0 06 

0 t 

0.15965015 10“^ 


0 2 

0 15905035 X 10"^ 


0 5 

0 15805068 X 10"“^ 


04 

0 15665114 X lO"’’ 



5 3 


He 


0 02 

0 19062364 

X 10~ 

0 04 

0 18221351 

X 

0 06 

0 17496617 

X ICT 

0 08 

0 16837958 

X ICT 

0 1 

0 16218441 

X 10” 

0 2 

0 13377557 

X kt 

0 3 

0 10668601 

X icT 

0 4 

0*79850001 

X icT 

0 5 

0.53748746 

X ICT' 

0 6 

0 30072266 

X 10”‘ 



Sable 5 4 


a 



o o 

0 1 

0 49270213 X icf^ 

0 2 

0 49055965 x icT^ 

0 3 

0 48S98883 x tCT^ 

0 4 

0 4Sl9®69 ^ ICT^ 

0 1 

0 47236286 x io"^ 

C 2 

0 46941452 x kT^ 

0 3 

0 46450CS9 X ICT^ 

0 4 

0 45762112 >' KT^ 

0 1 

0 45835926 X icT^ 

0 2 

0 45485608 x icT^ 

0 3 

0 44901745 x 

0 4 

0 44C«4356 X icf^ 


0.05 



WMmtOLAimE 


Corffieotiratioja of the solute in the fluid 
molecular diffusion coefficient 
equiralent dispersion coefficient 

hei^t of the rectangilar channel cnrer the central line 
reaction rate constant 

characteristic length along the axis of the channel or tube 
pre sure 

flux of the solute across a section of the clmnnel or ttibe 
radius of the tube 
radius of the core 

Telocity in the non-core region of the rectangular chanr»l 
velocity of "fee core in the reotacgular oham^ 

Telocity of the fluid relative to a plane movitB with the 
mean spe^ of the fluid 

rnmn speed of the fluid in rectangular channel 

velocity of the fluid in the mn-coare region of circular 
■fehe 

vaLociiy of the oolite in the flow ■&rou#i a circular tube 

velocity of the fluid in the circular tube, relative to a 
plane moving vd.th the mean speed of the flow 

mean speed of the flow in the circular tube 

hei^t of the core from the central line la the case 
of flow in a rectangular cbanael 



vxscosity of the fluid 

rBn-diaensiomlised. ayial distuuce relatiYe to a plane 
moving with the mean speed of the flow 

non- dlmensionalisi^ transverse distance 

non-dhaensioDBlised height of the core frc© Hie central 
line of the rectangular channel 

shear stress 


yield stress 
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CiiAPz:^ /I 


r m CGiTIEOOS EE^CTIi Oh VB LISEEESIC L i.A.’UllH t^(J/ 

or A BnrHLi fluie 

6 1 I TRmCTIOH 

xn. Cbapter IV ve have investigated tbe effects of viscosity 
variation on the dispersicm of a solute in the nresence jf hoiaogeneous 
ar»i heterogeneous reactions, in a fluid of variable viscosity, 
flowing between parallel plates and in circular tubes, using faylor' a 
approach In Chapter V, we have studied the dispersion in a Caason 
nod el fluid fLcwing in the two geometries stated above, using Baylor's 
laethod The effects of hoBc^neous reaction have also been investigated 

In this Chapter, we study a^an the effects of howagenecus 
reaction on dispersion of a solute in a jwn-lertonian fluid flowing 
betreen parallel plates aat^ in circular tubes by considemng the 
linghaia atodel 

6^ ET-— c~’ cr 1 "' ogtib''tj 3 if cttg? cr 31371715117 ii sr pl ctic 

L IE •ba li G rErVEEL' FARALLil F ^TES 

Gcwsider the lamnar flew of BinghaB fluid uisier a ccKjstant 
pressure gradient b^lween two parallel plates distant 2h apart ®ie 
physical configuration and the coordinate system are shewn in 
Fig (6 1 ) 

The ©luation of motion for one dioensional steady laminar 
flow of Bingtoa fluid is given by 




mmmi 


' ' ’ ' y ■ ,y yi 1 -' / ""'V i 

"/' '::'''' -y^y/yy''? f^’n 




';v' /■"^-^VviV?;;:;''-; 

WlliM& 








aM the constitutive equation is given bj 


(T>T^j (6 2-a) 

§ « 0 2-b) 

wtere Tq is the yield stress of the fluid and V is the ;-iscosity 
Integrating equation (6 1) and using the houndarj coiadition 


T =a 0 at y *» 0 


we get 




(6 5) 


If is the height of the core, we have 


T • y 


(6«4) 


Si i.Mtm ting T from ©qoations (6 2-a) and (6 5) sol using no slip 
condition at the wall after integratim, we obtain 


v- , 1 [f'^) g ' ^" ') - TgCh-^)! 


(6.5) 


in -^he region (y^ 1, y i. 


The velocity of the core is given hy 

(»-yo)5 
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ir t e region. (0^ j /^) 

■'■a'bstxtu^ing ti = j/h, rt^ = f^i' i” coucitxors (6 6) and (6 5) ani 
di/xdirg we obtain 

/ (n - 

7“ “ 1 ^2 (6 7-a) 

o (l - rig) 

(0<n<n_) (6 7-15) 

0 

The aveiTige velocxty Y which xs defined by 
h 

f = ^ / T dy (6 8) 

-h 

xs ol5t allied as 


"V* *» 




irhe 3 ?e 




(6 9 ) 

( 6 . 10 ) 


We now asstme that the dlffueing mtter while dispersing undergoes an 
inwersihle cteiaical reaxstion in. the flmd under isotbeisal ccraiitxons 
fhe diffusion ^uaticm in this case can 1 » written as 

2 

^ 'ix ** ^ "■ 

3 y 

where D* is the Kolecular diffusion coefficient and K is the reaction 
rate constant Following Taylor (l 953 ) aa the last two chapters, 
the differential equation governing the concentration of the solute 
relati-^ to a plane jsovxng with the oecn speed of the flow is giren hy 
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3^0 2 ^ 


(6 12 ) 


where 


2 Fh 
a * 


(6 15 ) 


aid. ¥ , the Telocitj of the fluid rt.lative to tte aovxng plane ib 
as follows 


- 1 

V „ V [(1 - 1) 2_, 

X '‘in / . 

m(1 - hg) 


¥ 4 - 1 ) 

'■a * 


(n < n < 1) 

O — — 


(0 ^ n £ Pq) 


(6 14-a'' 

(6 14-h) 


Ehe "boondaiy corditicais are 


(i) 

(ii) 


3c 

3o 

3n 


0 at n = 0 

0 at n e 1 
3 o 


(xii) o aM — are oontxmews at n » 


(6 15-e' 
(6 15 - 1 >) 
(6 15 - o ) 


Sol-vlng the equation (6 12) with the shore hoursiary ocMdxticms, we 
obtain tbe oonoentratiEm of tte solute aa follows 


o • H [ 1.^ ooEto an + B.J sinh an+itn)l (n^inll) 


where 


^ g 

®a (1-^) sxhh o 


etnh OTu cosho 

[(1-0 + r i 


Z sihh on, 


^1 “ TZ Y 
jaa (l-nQ) 


(6 16 ) 


(6 17 “^ ) 


(6 17-h) 
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F{n j 


[■ 


(n . 




- 2 /. ,2 

“ (1 - %) a (1 “ T3^) 


-"o)j 


H . 4 ^ (#) 


c » H [ip cosh an 5 1 f. n £, i^q) 

5 aa 


There 


^2 




■[ 


sinh (ot - an^) 


sinh <* 


C-' - 


(6.18) 

(6 19) 

(6 20 ) 


(6 21 ) 


when tIq = 0» w® case of a lewtonian fluxd for which c is 

deduced froa equation (6 l6) as 


o » oosh an 

U PL a sitii a 




(6 22 ) 


■srhioh has heen ohtsdmd by Gupta et al (1972) 

The TOluaetrio flew rate at which Ihe solute Is taaasjOTted 
aesroBS a seotiem. of ilis channel of unit hre^th is gxven by 


n 1 

0 « 2 J c Y dn + 2 / c Y dn 

^ X. X 


( 6 . 23 ) 


Substituting anprojaciate expressicsas for o and Y^ froia eqTsations (6 l6), 
(6 20), (6 t4-a), and (6 H-b) and integrating, tre get, 
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Q « (1°) . -y ... y T'(n ,a) 

m ot 


(6 24 J 


wnere 




( 1-1)3 (1-T\j)" 


A . 




„ ^T_ 8 Binh (®n ) 

4 coth g ^ 


o^Cl-rijj)^ o^Cl-no)"* sinii a 


4 sinh (on^) sxnh (a - an^^) 

______ _____________ 

a ( 1 - Tii_} sxiih a 


(6^5) 


Goiapiarxng with Picsk»s law of diffusxm we get toe ©<iuivalent diaperaion 
ooeffxcient D* as 


m 


t 


1 


B 2 2 
B a 


y(%»a) 


(6 26 ) 


'ib 0 n. n *• 0, we get B* as 
0 


D* SiL 

rj^O D 


^ 1 r1 1 . ooth g 1 1 

** '' ~ 'lr^ ^ ^ *" ""TT J 


^ 45 


2 '5' ‘ 


(6 27) 


5d a 


whxoh coittcides wxto the xesalt ohtainM hj Gupta et al (l972) 


^hen <» *»■ 0, i e toex® is no hc®ogenecw.s reaoticw., ijedTJoes to 




^ 16 ^ ''2+n„ ^ 


a*o B 105 -- -o 

which has heen ohtaiivsd hy Pan and Wgng (1966) 


(6 28 } 
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To stud/ the effe^^ts of dcmogeneous reaction on tne dispersion 
coefficient, P(n^,a) has been caleiilated and tabiilated in table (6 1 ) 

Table 6 1 


% 

a 



2 0 

0 13492921 X 10"’’ 


4 0 

0 73129225 x 10"^ 

0 01 

6 0 

0 41691002 X 10"^ 


8 0 

0 26109522 x 10"^ 


10 0 

0 17664779 x lO"^ 


2 0 

0 12957167 X lO"’’ 


4 0 

0.7C945273 X 10"^ 

0 05 

6 0 

0 40168828 x lO"^ 


8.0 

0 25186767 x icT^ 


10 0 

0 17054975 ICT® 


2 0 

0 12229157 X Id"’’ 


4.0 

0.66572481 X icr^ 

0 1 

6.0 

0.581 11 9C® X icT^ 


8.0 

0.25942817 x icT’^ 


10 0 

0.16254391 10"^ 


2 0 

0.10616444 X 10"*^ 


4 0 

0 3e239Sm X 10”® 

0 2 

6 0 

0 35584326 x 10*“^ 


8.0 

0.21212669 X 10"® 


10 0 

0.14437061 X iO“® 
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a 

--Cn„,a) 

o 

2 0 

0 88650938 X l<r^ 


4 0 

0 491875® x ICT^ 

0.3 

6 0 

0 2®67922 ICT^ 


8 0 

0 18251791 10"^ 


10 0 

0 12491804 X 10”® 


2 0 

0 70614697 ^ 10”^ 


4 0 

0 39797627 >« ICT^ 

0 4 

6 0 

0 23537860 x 10”^ 


8 0 

0 15151905 X icr^ 


10 0 

0 1(H51067 10’"^ 


2 0 

0.52917812 ICT^ 


4 0 

0.304404® ^ ICT^ 

0 5 

6 0 

0 t®55776 i(f^ 


8 0 

0 11995229 X 10*^ 


to 0 

0 83633334 * 10*"^ 


It am be seen fiom me above table that for a givaa n^iHn^,a) 
decreaaes m a iocsreaees and it decreases further as 

6 3 EEODI OF HCWOiSEOUS REACTION OR OISPEESIDR IR IMIHAE 

oj bbigelam; hbib i j a ciecdlar rosr 

In this section we study liie effects of ho®ogeneoua reaction 
on dispersiou of a solute in the laminar flow of a Bin#aa fluid in a 
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cixonlex ti te The paY'-xCdl omf’iguraticm ard coordinate system 
are siorm in Fig (6 2) The equation of irntiun for one dimensional 
steady laminar flo7 of Bingham fluic is gxren tj 


dx r dr ^ ^ 


(6 29) 


p,rtr! the constitutive equation for the Binglais fluia is 

K-g)-t-to (6 30^1) 


dr 


0 


Ct 5 


(6 50 ~h) 


where is the yield stress and v is the viscosity Following 

toB slailax procedure as in tl» last seotxwi, the velocity profile 
is given hy 

V . -i K-H) (^) - T„(E-r)] (6 31) 


in the region(rQ ^ i. ^) » 


B 


t(-S) 


(6 52 ) 


m the xe^ion (O ^ r £. 3 ^ 0)1 'whai^ the radius of the oorsj is 
given hy 


2^ 

(-dp/dxy 


(6 35 ) 


Juhstittttlng p » r/^ and 


profile an 


T « T [1 - 
o 


(n - n 

2^1 

(1 - 


^0 » we can write the velocity 

(*>0 1 1 1 't) 54-a) 

(0 < n < n ) (6 34-12) 

V — T— O' 






■ : ' 

5 ' ^1) ^ t- ; 


-“■ .' 5 !’^{ ? 

sf 


?■■: . '- 
■ii^iiii 
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Th 0 arerage Telocit/ / j-S as 

_ H 

^ / r ^ dr 

R 0 

Using tl^e equations (6 34“S-)> (6 34 ""b) m equa'*"-.c*i 
aTerage Telocity as 


(6 55) 

(6 55) ge'b 


Y = 




(6 56) 


where 


n 


5 + 2ng + 

g 


(6 57) 


•Phe velocity of the fluid relative to a plane Boving with the mean 
apeed of the flow is obt^ned as 


(n - n ) 

** ri-n ^ ^ 1 

Y » T - Y » ¥ [-^ - — - :tJ 

n(l -Hq) 


(6.58) 


in the region ( 5. 1 ) > and 

Y ja ^ 59} 

X h 

in, the region ( 0 

As hefore, the diffusion equaticaa in this c&m is gtrm hy 


aJL 

r 3x 




at 


Q » 




m+ 1£ 

ax^ St 


(6 40 ) 


rollcmrjg Taylor (1953) the non-dismsionalie^ Aiffiiatcn e^uatim 
goreming the conoertratxcn of the solute rela-i/e to a plane moving 
wit-0 libe lasaji apssd of tiio fleas' is obtain^ as 
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ac 


(n - a c » H [0 - 1) - 




n ar) 3 rj^ ----- ivj, - -/ - - _ in 1 1 ) 


n(i - n J 


1l-r) 


(Oln < n ) 


(6 4l-a) 
(6 41-'b) 


"here 


5 » S— (^) Y 
“ DL ^ 35 ' 

The boujndary cotri-itiom are 

^i) |~ = 

on 


an 

c aM •“ 
an 


at n = 0 

(6 42 -a) 

at n = 1 

(642-h) 

are contimious at n ® n„ 

o 

(6 *42 -o) 


Solving the eqtiatiO!t»(6*41*^) awl (6 41 "l^) with the above hoancLary 
ooiiiiti<ms, we ohtaxn itse ooncentratioa profile as follows 


o - (an) ^ SLLf^ (0 < n in^) 


(6»43) 


n a 


where 


h(l - ng ) 


H r/^\ 4 > ^ 1 ^“^ I^(ano)K^(a) 

— 3 - i(_+3-}n t rrri ^ 


x;(^ 


*_ ( — ^ 


I^(a) 


( a) 


( 644 ) 


«/ II^(<*)K^{at) + I^(ot) K^(a)l dt 


(6 » 45 ) 
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aM dre :^odxf_ed ^sc-1 far:c+xona of f.rst and cecond kind of 

Older zero ard are of order 1, and 


c » Ag ijoin) + Kjan) - 


H 


[ 




^(1 - n^) 


C, 

-(%+%-2n^ F(n)] (n <n<i) 

or <x 


(646) 


wliere 




(6 47) 


:'(n) » / [IqCw) K^(at) - I^(at) K^Con)] ^ dt 


(6 # 48 ) 




-H 


n(l -n^) 


2n„ . n„ I.(cm„) 

2 [KiC») - 


y XAa) 


2 T{ 1 (m ) 

"o <>■',, +^-2_ + 


a I,(o) 




0,5 


( 649 ) 


-H n. 2C , 2n I (w ) 

^ ^ M o o o 


\ ^ [(—+%) I^(an^) - 

nCl - %) 


(6 50) 


iTq ** Of we get ttie osae of a l@irtom.aa fluid for wbioh we 
obtain Of from equation (6 46) by putting « 0, as 


-4H lj«n) 


, 2H n , 9k _ I_ 


(6 51 ) 


•jidiich is the saffls as obt aan^ earliCT* 
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35'o TOliffiietric irate at iriaica tae solute is tran^jortea across 
a section of liie tube is given b? 


.2 A 


1 


w = 2wR f ° n c dn+2wK J n c 


(6 52 


o 


aibavxtutxHg appropriate axpresaions for c and frcm 
eqtuatxons (6 45)f (6 46 )> (6 3^) and (6 39) xn ihe eqmtxons (6 52) 
we obtain 


Q_ _ 26^ fl£\ r ^1 

2 ^ Ti ^ o 


vB 




D 3X'' 


2n^ 


] 


(6 53^ 


ifeere 




2 A^(i-.n) n “ Hq 

Hi 


(6,54) 




{Hq I^CoiIq) - I^{a)) + { 


I^(a) aigCa) 


a 


a 


iXon ) 2n^ ^ ®2 r / ^ 

_2 — J — 2^ + Q— I — 2-) +>^ (K^Ca) 

1 (a) mJa) 
Tj^EXan )> - e_ {— — + 5— 


-2 ‘ a 


T? K (an_) 2n^ E (cm ) (l - nb 

o I o o d. o ^ , y_ 

5 i 2 

^ 2a 


«(l . n*) 2« Cl - T,®) 2n„ *(l - n’) 

11^ III ^ ' 


2a 


4 


3a 
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5a^ 3 q ? 6 a ^ / 

1 2 2 
- 2t1q £^ / Ti lo(an) an - 2ng e^ J n ’"gCan) dn 


• ’Z 

- 2n / (^n + 2n_ n -• n ) ^(n) dn 

O ^ W 

and 

^2 ^ (l — tIq ) 

^1 i 

Bg n (l - T ^)^ 

^2 ” i 


(6 55 ) 


(6 56 ) 


(6 57 ) 


A., A , B f d » ^(n) axe given iw eqmtiona <6 44), (6 49), (6 50), 

I 2 2 

(6*47) aM (6,48) 


Garnpaxlm «ifeb Ixck’s law of dtffasion, to get tbe ajtiiTalerrt dieperaion 
oo^Xlcieat aa 


*a 


2 n o " n^(l - nn ^ 


2 „2 


(6 58 ) 


Iben Tj^ 3= 0, TO get Ifae caae of a ^wtonlsa fluid focr tiiiicfa 


JP la gijreft tj 

ill M g 

[-c + — o jl 

^=0 B 5- / \ 4 

o I ^ Va > 3 a a 


(6 59 ) 


which XB liie same as obtained to C&apter (iv), ecixatlon (4.60 ) 


SIq stiud^ 14ie effects of tbe hosaos^csous refsctioa oct the 
equivalent aa£^«rsxon coefficient, equation (6 58) ie slBqiLxfted for 
t w!fn\X n^f n^ectlr® and tiller powers, as 
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m - 


2P 




o 4, [C c 5 

D n (1 - Hq) a I^(o) 24a 


412(a) 


a4 


{ 

*■ ‘*0 a 


S 4S I„(a) 

h 1 2 1 

a ”1(5 


l5o 


_8 

3a^ 


T ^ ^ 


2855 


a (a ) 




a2 


5 20 448 20736 


■>}J 


zrH 

F 


ir.2 


^(n^,a) (6 60) 


It xs oDoervei 1iiat as a 0, eqxxatioa (6 60 ) reduces to 




5 44n 

r 2 o-i 

I A - 35 i 


which has been obt aine d by Ian and Wang (l966) 


(6 61 ) 


Bjiiatian. (6 60} tea been calculated and tabulated in table (6 2/ 
£jxm #i.ioh it can be seen that lixe equivalent d.i£5>«rslon coefficient 
dBoieases as a Increases and it decreases further ’utoen al«> 


increases 
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l?3fele6 2 


0 00 


0 04 


0 08 


0*1 


H&re 


o 

0 04 
0 1 
0 3 
0 5 

0 04 

0 1 
0 3 
0 5 

0 04 
0 1 
0 3 

0 5 

0 04 
0 1 
0 3 
0 5 

F(T)^,a) is the coefficient of 


a) 

0 2CB51111 10""® 

0 208l94‘t4 X icf"® 
0 20708333 10"^ 

0 204^111 X iCf"* 

0 20799367 x 10"^ 
0 20787716 »< ICT"' 
0 20676746 x icf"® 
0 20454a» ’< 10"^ 

0 20767624 x ICT^ 
0 20755987 ICT"* 
0 ac«45t59 X ICf^ 
0 204235<» X tcT^ 

0 OT51753 X id""® 
0 20740123 X icf"® 
0 20629365 x 
0 20407348 x lO""® 

?.n2 

( — g — ) in equation (6 60 ) 



164 


6 4 (X)2fgi;^CIOI<3 

SSie effect of htMogeiKous reaction on (ii^ersion of a aslute 
in laainar flow of a Bin^aa fliiid xn rectargular cbeinnel and circular 
pipe has been studied by using $^lor»s i^proach It ia ob ©erred liiat 
tae equivalent dispersion coefficient decreaBes as the reaction rate 
constant mcreasesj for a given However it d^c^aaea furtiier if 
increases 

It may be reoafteed here that the non-Hewtonian beiiaviaur of 
the fluid decreases the Baylor's diffusion coefficient for a given mean 

velocity. 
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c cone eafcrat ion of a aslute in ihe flowir^ solrent 

B aolecular oiffusion coefficient 

a*- equivalent (liffusion coefficient 

h hei^ t of tiie rectangular ctonnel froa the central lone 

K reactior rate constant 

L chaxacterxBtic ler^h along the tube or charnel 

p pressure 

Q flux of t^e ailute across a section of the cfcannel or tube 

1 radius of liie circailar tube 

r radms of 14ie core in the case of flow in circular tube 
o 

¥ Telocity in the non-core region 

¥^ Telocity of -Sie core 

y^ bei#it of the core traa tne central line 
p Tisooelty of ttie fluid 

5 nim-gioenaiomlis^ axial length relatlTe to the both® plane 

n nonudiaeasionallsed ler^h in the transrsrae direotion 

T sh®ir stress 

V yield stress of the fluid 
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